US  Army  Corps 
of  Enginears 

Waterways  Experiment 
Station 


Technical  Report  HL-96-7 
August  1996 


96-02178 

llllllllll 


i 


Return  Velocity  and  Drawdown 
in  Navigable  Waterways 


by  Steve  Maynord 


AD-A286  906 

an  iM  nni  ntl  IlH  HI 


Approved  For  Public  Reloase;  Distribution  Is  Unlimited 


Prepared  for  Headquarters,  U.S.  Army  v^orps  of  Engineers 


The  contents  of  this  report  are  not  to  be  used  for  advertising, 
public, ation,  or  promotional  purposes.  Citation  of  trade  names 
does  not  constitute  an  official  endorsement  or  i^roval  of  the  use 
of  such  commercial  products. 


PRINTED  ON  RECYCLED  PAPER 


Technical  Report  HL-96-7 
August  1996 


Return  Velocity  and  Drawdown 
in  Navigable  Waterways 

by  Steve  Mayrrard 

U.S.  Army  Corps  of  Engineers 
Watenways  Experiment  Station 
39(^  Halls  Ferry  Read 
Vicksburg,  MS  39180-6199 


Final  report 

Approved  for  public  release;  distribution  is  unlimited 


Prepared  for  U.S.  Army  Corps  of  Engineers 
Washington,  DC  20314-1000 


Waterways  Exparimant  Station  Cataloging’in-Publlcatlon  Data 
Maynord,  Stephen  T. 

Return  velocity  and  drawdown  in  navigable  waterways  /  by  Steve  Maynord ; 
prepared  for  U.S.  Army  Corps  of  Engineers. 

102  p. ;  ill. ;  28  cm.  —  (Technical  report ;  HL-96-7) 

Includes  bibliographical  references. 

1 .  Channels  (Hydraulic  engineering)  —  Mississippi  River.  2.  Barges  —  Missis¬ 
sippi  River.  3.  Navigation  —  Mississippi  River.  4.  Drainage  —  Mississippi  River. 
I.  United  States.  Army.  Corps  of  Engineers.  II.  U.S.  Army  Engineer  Watenways 
Experiment  Station.  III.  Hydraulics  Laboratory  (U.S.  Army  Engineer  Watenways 
Experiment  Station)  IV.  Title.  V.  Series:  Technical  report  (U.S.  Army  Engineer 
Waterways  Experiment  Station) ;  HL-96-7. 

TA7W34no.HL-96-7 


Contents 


Preface  . vii 

1—  Introduction  .  1 

Background .  1 

Objective  and  Scope  .  2 

2—  Application  of  Confined  Channel  Methods  to  Large  Waterways  .  3 

Previous  Confined  Channel  Studies .  3 

Modification  of  Confined  Channel  Methods  to  Large  Waterways .  5 

Large  Waterway  Studies .  5 

3—  Return  Velocity  and  Drawdown  Data  Sources,  Errors,  and  Analysis  .  .  9 

4— 1991  and  1995  Analytical  Methods . 12 

Return  Current  Comparisons  with  1991  Analytical  Method  . 12 

1995  Analytical  Method . 13 

5—  Development  of  1995  Analytical  Method  for  Prediction  of 

Water-Level  Drawdown  . 47 

6—  Comparison  of  1995  Analytical  Method  with  Independent  Field  Data  .  .  52 

7—  Program  NAVEFF . 72 

8—  Results  and  Conclusions . 86 

References  . 88 

Appendix  A;  Notation .  A1 


List  of  Figures 


Figure  1.  Navigation  effects;  A  =  return  velocity;  B  =  drawdown;  C  = 
propeller  jet;  D  =  wake  flow;  E  =  bow  wave;  and  F  =  slope 
supply  flow .  1 

Figure  2.  Cross  section  of  the  Illinois  River  at  the  Kampsville  site  for 

tripl.  . 11 

Figure  3.  Observed  physical  model  return  velocity  versus  computed  return 
velocity  using  1991  analytical  method . 13 

Figure  4.  Observed  prototype  return  velocity  versus  computed  return 

velocity  using  1991  analytical  method . 14 

Figure  5.  Observed  physical  model  average  return  velocity  versus 

computed  Schijf  return  velocity . 16 

Figure  6.  versus  using  physical  model  data . 17 

Figure  7.  V^/V„  versus  N,iie  usipg  physical  model  data . 18 

Figure  8.  Observed  Kampsville  physical  model  return  velocity  versus 

computed  return  velocity  using  :i>55  analytical  method . 20 

Figure  9.  Observed  Kampsville  physical  model  return  velocity  versus 
computed  return  velocity  using  1995  analytical  method  for 
individual  tests . 21 

Figure  10.  Observed  Kampsville  physical  model  return  velocity  versus 
computed  return  velocity  using  1995  analytical  method  for 
individual  tests  . 23 

Figure  1 1 .  Observed  Kampsville  physical  model  return  velocity  versus 
computed  return  velocity  using  1995  analytical  method  for 
individual  tests  . 25 

Figure  12.  Observed  Kampsville  physical  model  return  velocity  versus 
computed  return  velocity  using  1995  analytical  method  for 
individual  tests  .  27 

Figure  13.  Observed  Kampsville  physical  model  return  velocity  versus 
computed  return  velocity  using  1995  analytical  method  for 
individual  tests  . 


29 


Figure  14.  Observed  Clark's  Ferry  physical  model  return  velocity  ve-sus 
computed  return  velocity  using  1995  analytical  method  for 
N<52 . 33 

Figure  15.  Observed  Clark’s  Ferry  physicpJ  n«)dt  return  velocity  versus 
computed  return  velocity  using  l'-";  .  oalytical  method  for 
N  =  85 .  34 

Figure  16.  Observed  Kampsville  protOi, ; velocity  versus  computed 
return  velocity  using  195:j  OiiSiiytiea'  method . 35 

Figure  17.  Observed  Clark’s  Ferry  ,  .ot  xype  return  velocity  versus 

computed  return  velw  i  ly  ttsing  1995  analytical  method . 36 

Figure  18.  Observed  Kan^sville  pi  otti* /pc  return  velocity  versus  computed 
return  velocity  using  1!^>95  k  u  l^^ical  method  for  individual  tests  37 

Figure  19.  Observed  Kampsville  prt'to;  .pc  return  velocity  versus  computed 
return  velocity  using  1993  analytical  method  for  individual  tests  39 

Figure  20.  Observed  Kampsville  prototype  return  velocity  versus  computed 
return  velocity  using  1995  ana’.ytical  method  for  individual  tests41 

Figure  2 1 .  Observed  Kampsville  prototype  return  velocity  versus  computed 
return  velocity  using  1995  analytical  method  for  individual  tests43 

Figure  22.  Observed  Kampsville  physical  model  drawdown  versus  computed 


drawdown  using  1995  analytical  method . 49 

Figure  23.  Observed  Clark’s  Ferry  physical  model  drawdown  versus 

computed  drawdown  using  1995  analytical  method  . 50 

Figure  24.  Observed  Kampsville  prototype  drawdown  versus  computed 

drawdown  using  1995  analytical  method . 51 

Figure  25.  Observed  ESE  prototype  return  velocity  versus  computed  return 
velocity  using  1995  analytical  method  . 53 

Figure  26.  Cross  section  for  ORL  tests  89192-89196  .  54 

Figure  27.  Cross  section  for  ORL  tests  89198-89203  .  55 

Figure  28.  Cross  section  for  ORL  tests  90172-90179  .  56 

Figure  29.  Cross  section  for  ORL  tests  90191-90195  .  57 


Figure  30.  Observed  ORL  prototype  return  velocity  versus  computed  return 
velocity  using  1995  analytical  method  . 70 


Figure  31.  Observed  ORL  prototype  drawdown  versus  computed  drawdown 


using  1995  analytical  method . 71 

Figure  32.  Listing  for  PC  program  NAVEFF . 73 

Figure  33.  Output  from  PC  program  NAVEFF  . 85 


Tables 


Table  1.  Physical  Model  Experimental  Conditions  . 31 

Table  2.  Kampsville  Prototype  Experimental  Conditions . 45 

Table  3.  Observed  Versus  Computed  Drawdown  for  WC  Norman  Test  51 

Table  4.  ORL  Test  89192-891% . 58 

Tables.  ORL  Test  89198-89203  .  61 

Table  6.  ORL  Test  90172-90179  .  64 

Table  7.  ORL  Test  90191-90195  .  66 


vi 


Preface 


The  work  described  in  this  rq>ort  was  conducted  at  the  U.S.  Army 
Engineer  Waterways  Experiment  Station  (WES)  during  the  period  June  1995 
to  September  1995.  This  investigation  was  sponsored  by  Headquarters,  U.S. 
Army  Corps  of  Engineers  (HQUSACE),  under  the  Navigation  Hydraulics 
Research  Program  as  part  of  the  Civil  Works  Investigation  Work  Unit, 

"Vessel  Generated  Forces  and  Protection,"  under  HQUSACE  Program 
Monitor  Mr.  Sam  Powell. 

The  study  was  accomplished  in  the  Hydraulics  Laboratory  (HL)  under  the 
direction  of  Mr.  Richard  A.  Sager,  Acting  Director;  Mr.  Robert  Athow, 
Acting  Assistant  Director;  and  Dr.  Larry  L.  Daggett,  Acting  Chief,  Naviga¬ 
tion  Division  (HN).  The  study  was  conducted  and  the  report  was  written  by 
Dr.  S,  T.  Maynord,  Navigation  Effects  Group,  HN. 

At  the  time  of  publication  of  this  report,  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commande,  was  COL  Bruce  K.  Howard,  EN. 


Thr  contints  of  this  report  are  to  be  used  for  ad^rtuutg,  publication,  or 
promotional  purposes.  Otuticn  of  trade  names  does  not  constitiire  an  official 
endorsement  or  approval  of  the  use  of  such  commercial  produas. 


1  Introduction 


Background 

A  vessel  moving  in  a  navigation  channel  displaces  water  and  creates  a  flow 
of  water  (return  velocity)  alongside  the  vessel  in  a  direction  opposite  that  of 
the  vessel.  The  magninide  of  the  return  velocity  is  primarily  dependent  on 
vessel  speed,  channel  cross  section,  and  the  submerged  cross-sectional  area  of 
the  vessel.  The  return  velocity  is  also  accompanied  by  a  lowering  of  the 
water  level  alongside  the  vessel,  drawdown,  as  well  as  other  navigation 
effects  shown  in  Figure  1.  Return  velocity  and  drawdown  have  been  exten¬ 
sively  studied  in  prismatic  cross  sections  and  confined  or  restricted  channels, 
i.e.,  those  having  a  low  blockage  ratio  N,  the  ratio  of  channel  cross-sectional 
area  to  vessel  cross-sectional  area.  Previous  studies  (l^rmanent  International 
Association  of  Navigation  Congresses  (PIANC)  1987,  Schijf  1949)  have 
addressed  the  channel  subility  aspects  of  return  velocity  and  drawdown  in 


Figure  1 .  Navigation  effects:  A  =  return  velocity;  B  =  drawdown;  C  •= 
propellor  jet;  D  »  wake  flow;  E  =  bow  wave;  and  F  »  slope 
supply  flow 
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confined  channels.  In  other  studies  (Blaauw  and  van  der  Knaap  1983,  Tothill 
1966),  a  primary  concern  of  drawdown  has  been  the  resulting  lowering  or 
squat  that  the  vessel  experiences  when  drawdown  occurs.  In  confined  chan¬ 
nels,  vessel  speeds  can  be  limited  by  vessel  squat  when  depth/draft  ratios  are 
low. 

In  a  general  sense,  return  velocity  and  drawdown  are  analogous  to  pier  or 
obstruction  effects  in  a  river.  For  piers  that  are  small  relative  to  the  size  of 
the  river,  large  effects  are  found  near  the  pier  and  negligible  effects  near  the 
bank.  If  the  pier  took  up  a  greater  portion  of  the  river,  the  effects  would  be 
large  near  the  pier  and  significant  near  the  bank.  Finally,  as  the  pier  begins 
to  dominate  the  cross  section,  the  effects  near  the  pier  and  near  the  bank  are 
similar  in  magnitude.  Because  return  velocity  and  drawdown  follow  these 
trends,  mathematical  decay  functions  whose  shape  varies  with  the  vessel  cross- 
sectional  area  relative  to  the  cross-sectional  area  of  the  river  are  used  to 
describe  the  distribution  of  return  velocity  and  drawdown  between  the  vessel 
and  the  bank  line.  For  a  large  vessel  size  relative  to  the  river,  the  decay 
functions  should  collapse  to  an  almost  uniform  distribution  from  bank  to 
vessel. 

Return  velocity  and  drawdown  are  of  interest  in  assessing  the  effects  of 
navigation  on  environmental  concerns  for  large  navigable  waterways. 

Channel  stability  and  vessel  squat  impacts  due  to  return  velocity  and  draw¬ 
down  are  often  small  in  large  waterways  like  the  Ohio  and  Mississippi  where 
blockage  ratios  become  large.  One  exception  is  areas  where  vessels  travel 
near  bank  lines  such  as  in  lock  approaches.  The  analytical  techniques  for 
estimating  return  velocity  and  drawdown  developed  for  confined  channels  are 
not  directly  applicable  to  large  waterways  because  many  basic  assumptions 
used  in  their  development  are  not  met  in  large  waterways. 


Objective  and  Scope 

This  study  develops  techniques  for  estimating  return  velocity  and  draw¬ 
down  from  commercial  vessels  operating  in  navigable  waterways  and  com¬ 
pares  them  to  model  and  prototype  data.  Confined  channel  methods  are 
modified  to  extend  their  application  to  large  riverine  navigation  systems.  The 
results  are  incorporated  into  the  PC  program  NAVEFF  with  this  report 
serving  as  the  documentation.  This  study  does  not  address  other  navigation 
effects  such  as  waves  or  propeller  flows  shown  in  Figure  I .  This  report 
presents  a  modification  a^  an  expansion  of  the  results  presented  in  Maynord 
and  Siemsen  (1991)  and  uses  data  from  shallow  draft  navigation  on  the  Ohio, 
Illinois,  and  Upper  Mississippi  Rivers. 
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2  Application  of  Confined 
Channel  Methods  to  Large 
Waterways 


Previous  Confined  Channel  Studies 

An  excellent  review  of  techniques  to  detemtine  squat,  drawdown,  and 
return  velocity  in  confined  channels  is  presented  by  Blaauw  and  Van  der 
Knaap  (1983).  Only  th*  rtost  pertinent  studies  are  discussed  in  the  following 
paragraphs. 

Schijfs  (1949)  conservation  of  energy  iqiproach  is  frequently  used  in  deter¬ 
mining  return  velocity  and  drawdown  in  confined  channels.  The  basic 
assumptions  required  in  this  approach  are:  uniform  trapezoidal  or  rectangular 
cross  section,  uniform  return  current  velocity  from  vessel  to  bank  line, 
uniform  drawdown  from  vessel  to  bank  line,  friction  losses  disregarded, 
uniform  or  negligible  ambient  velocity;  and  vessel  on  channel  center  line. 
Considering  all  motions  relative  to  the  vessel,  continuity  requires  that 


VA,  »  A,(V+V^' 


(1) 


where 

V  =  vessel  speed 
Af  =  undisturbed  channel  area 

==  disturbed  channel  area  around  midsection  of  vessel  excluding  A. 
V,  =  average  return  velocity 


'  For  convenience,  symbols  and  abbreviations  are  listed  in  the  notation  (Appendix  A). 
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=  submerged  cross-sectional  area  of  the  vessel  at  the  midsection 
Conservation  of  energy  requires  that 


^  (V  + 

—  +  h  =  — - - 

2g  2g 


+  (h  -  z) 


(2) 


where 

g  =  gravitational  acceleration 
h  =  undisturbed  average  water  depth 
z  =  average  water-level  drawdown 

The  unknown  quantities  V,  and  z  can  be  determined  from  Equations  1  and  2. 
Differentiation  of  Equations  1  and  2  leads  to  a  maximum  speed  (called  the 
critical  or  limiting  speed  Vj)  which  can  not  be  exceeded  by  a  self-propelled 
vessel  and  has  been  verified  in  both  model  and  prototype  investigations. 
PIANC  (1987)  presents  a  coefficient  0  that  varies  with  V/V^,  in  the  Schijf 
equations  to  improve  agreement  between  observed  and  computed  return 
velocity  and  drawdown.  The  energy  equation  is  rewritten  as 


2g"  ~  2g 


*  {h  -  z) 


(3) 


Gates  and  Herbich  (1977)  included  the  displacement  thickness  from 
boundary  layer  concepts  to  determine  effective  vessel  draft  and  beam  improv¬ 
ing  the  agreement  between  Schijf  equations  and  prototype  measurements  of 
vessel  squat.  Inclusion  of  boundary  layer  concepts  partially  addresses  the 
assumption  requiring  that  friction  losses  be  disregarded. 

Bouwmeester  et  al.  (1977)  reported  that  Schijf  (1949)  cannot  be  used  even 
approximately  for  N  greater  than  about  33  because  water-level  drawdown  and 
return  velocity  are  greatly  nonuniform.  Many  tows  in  large  waterways  have  a 
value  of  N  greater  than  33. 

PIANC  (1987)  documented  techniques  for  using  the  Schijf  (1949)  equations 
in  channel  stability  investigations  of  prismatic  channels.  Schijf  is  recom¬ 
mended  for  BJB  from  2  to  12,  where  is  the  channel  width  and  B  is  the 
beam  of  the  vessel.  Additionally,  the  ratio  of  channel  width  to  vessel  length  L 
can  be  a  significant  factor  in  defining  the  ratio  of  maximum  water-level  draw¬ 
down  to  average  water-level  drawdown  and  maximum  return  velocity  to  aver¬ 
age  return  velocity.  Channels  having  a  large  ratio  will  have  a 


Chapter  2  Application  of  Confined  Channel  Methods  to  Large  Waterways 


1 


well-defined  return  velocity  and  drawdown  time-history  during  the  time  the 
tow  is  aqjacent  to  a  given  point  on  the  bank. 
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Modification  of  Confined  Channel  Methods  to 
Large  Waterways 

Some  of  the  six  basic  assumptions  for  applying  the  Schijf  equations  are 
often  not  met  in  large  waterways.  The  ambient  and  return  velocities  and  the 
water-level  drawdown  are  not  uniform,  the  cross  section  is  not  prismatic,  and 
the  vessel  is  not  on  the  channel  center  line.  The  equations  determine  average 
values  of  return  velocity  and  drawdown.  Maximum  values  near  the  shore  are 
of  the  most  interest  for  both  channel  bank  stability  and  environmental 
concerns.  Environmental  interests  often  also  require  the  variation  in  return 
velocity  and  drawdown  between  vessel  and  shoreline.  The  techniques 
presented  here  predict,  at  any  given  point  between  vessel  and  shoreline,  the 
maximum  deviation  from  ambient  conditions  due  to  return  velocity  and  draw¬ 
down  during  the  tow  passage. 

A  major  problem  that  must  be  overcome  in  evaluating  asymmetric  channels 
found  in  large  waterways  like  the  Ohio,  Mississippi,  or  Illinois  Rivers  is  how 
to  handle  tows  off  the  channel  center  line.  One  technique  assumes  a  mirror 
image  on  each  side  of  the  tow  and  computes  average  return  velocity  and  draw¬ 
down.  This  is  satisfactory  until  the  tow  gets  near  one  bank  line  and  far  from 
the  other  resulting  in  a  substantial  amount  of  the  return  flow  passing  around 
the  front  of  and  beneath  the  tow.  Then,  the  return  current  predicted  on  the 
side  near  the  bank  using  the  mirror  image  channel  is  much  larger  than  actually 
occurs. 

A  second  problem  with  the  Schijf  equation  is  in  certain  asymmetric  channel 
shapes.  Any  channel  with  wide  shallow  areas  on  either  or  both  chatmel  sides 
is  particularly  subject  to  overestimating  the  return  velocity  at  high  vessel 
speeds.  This  problem  becomes  significant  when  tow  speeds  exceed  about 
90  percent  of  Schijf  for  the  total  section  and  when  the  ratio  of  maximum 
depth  to  average  depth  exceeds  about  1.3.  Fortunately  most  tows  are 
traveling  at  less  than  0.9  V^,  normally  from  50  to  75  percent  of  V],. 


Large  Waterways  Studies 

Hochstein  and  Adams  (1989)  documented  a  technique  for  estimating  return 
velocity  on  the  upper  Mississippi  River  that  is  used  by  Environmental  Science 
and  Engineering  (ESE)  (1981)  and  Simons  et  al  (1988)  and  on  the  Ohio  River 
by  U.S.  Army  Engineer  District,  Huntington  (1980).  The  average  return 
velocity  is  defined  as 


4 


I 


9 


Chapter  2  Application  of  Confined  Channel  Methods  to  Large  Waterways 


I 


V,  =  V[{aB^  -  flr  ♦  1)“  ’  -1] 


(4) 


4 


f 

I 


4 


1 
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where  a  =  (N/(N-1))*  * 

For  vessel  speeds  less  than  0.65Vt,,  where  Bj-  is  a  coeificieut  defined  as  Bf 
=  0.3  exp(1.8V/VL)  and  for  vessel  speeds  equal  to  or  greater  than  0.65Vl, 

Bf  =  1.  From  Maynord  and  Siemsen  (1991),  the  Schijf  and  Hochstein  aver¬ 
age  return  velocity  equations  give  similar  results  for  vessel  speeds  equal  to  SO 
to  60  percent  of  Schijf  Vi.  In  Equations  5  through  7,  Hochstein  presents  a 
method  for  determining  the  return  veiocit>’  distribution.  A  channel/tow  width 
factor  is  defined  as 


B 

a  =  max  (1,0.114  +  0.715) 

B 


(5) 


The  maximum  return  velocity  V„  at  the  vessel  is  determined  from 


=  a  V 


(6) 


For  a  -  1 ,  the  distribution  of  velocity  is  uniform  from  vessel  to  bank 
line.  For  1  <  a  <  1.5,  the  distribution  of  return  velocity  is  linear.  Foro; 

^  1 .5,  the  distribution  of  return  velocity  is  linear.  For  a  k ,  the  distribution 
of  return  velocity  from  vessel  to  bank  line  is  defined  as 


VJiY)  =  exp  i-Y/k) 


(7) 


where 


V^(Y)  =  maximum  return  velocity  during  tow  event  at  distance  Y  from 
center  line  of  vessel 

k  =  B.ifc/(a{l-exp[-F(a)al}) 

^lide  -  distance  from  vessel  center  line  to  shoreline  • 

F(a)  =  0.42  +  0.52  In  a 
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ESE(1981),  Simons  et  al.  (1988),  and  USAED  Huntington  (1980)  results 
showed  that  the  Hochstein  equations  (Equations  4  through  7)  underestimate 
measured  return  velocity. 

Maynord  and  Siemsen  (1991)  presented  a  method  for  estimating  return 
velocity  in  large  waterways.  Average  return  velocity  is  computed  using  Schijf 
(1949)  without  the  correction  factor.  To  handle  tows  navigating  off  the  chan¬ 
nel  center  line  in  asymmetric  channels,  the  average  return  velocity  is  then 
proportioned  on  each  side  of  the  vessel  using  the  following  equation  developed 
from  physical  model  studies  reported  in  Maynord  (1990) 


-^  =  C,  SKEW  *  q 


(8) 


where 

V„  =  average  return  velocity  for  each  side  of  vessel 
SKEW  =  A,/(2A,a,) 

=  cross-sectional  area  from  tow  center  line  to  bank  line 

Maynord  and  Siemsen  (1991)  used  coefficients  of  C,  =  0.36  and  (Zj  =  0.64 
in  Equation  8.  The  ratio  a’  of  the  maximum  return  velocity  on  each  side  of 
the  tow  V„„  to  V„  is  determined  from 


^  =  a'  =  max  [1,0.0247V^,  +  0.734J 


(9) 


where  =  2A,i^/A„. 

Note  the  similarity  between  Hochstein’s  Equations  5  and  6  which  use  a  width 
ratio  and  Equation  9  which  uses  an  area  ratio.  As  the  area  ratio  (N,;,),) 
increases,  the  distribution  of  return  velocity  becomes  nonuniform.  From 
Equation  9,  less  than  11  gives  a'  =  1  and  the  distribution  of  return 
velocity  between  vessel  and  shoreline  is  uniform.  The  return  velocity  distri¬ 
bution  from  tow  to  bank  line  is  determined  from 


B 


B. 


sid€ 


B 


)] 


(10) 


where 
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C=1.2(a’  -  1) 


(11) 


V^(Y)  is  the  maximum  return  velocity  during  the  tow  event  for  a  given  distance 
Y  from  the  vessel  center  line.  V,  (Y)  is  linearly  added  to  ambient  currents  for 
upbound  tows  and  subtracted  from  ambient  currents  for  downbound  tows.  Y 
must  be  greater  than  B,  because  from  Y  =  0  to  Y  =  B  the  velocity  is  also 
affected  by  flow  under  the  barges  which  can  be  considerably  larger  than  the 
return  velocity.  Based  on  a  limited  data  set  used  in  Maynord  and  Siemsen 
(1991),  the  Hochstein  method  underestimated  the  measured  return  velocities 
while  the  Maynord  and  Siemsen  (1991)  method  overestimated  the  measured 
return  velocities. 

An  obvious  limitation  of  these  distribution  equations  is  that  the  distribution 
does  not  vary  with  the  shape  or  local  depth  of  an  asymmetric  cross  section, 
although  both  methods  do  address  a  vessel  off  the  channel  center  line.  For 
the  remainder  of  this  report,  Equations  1,2,  and  8  through  11  are  referred  to 
as  the  1991  analytical  method. 
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3  Return  Velocity  and 

Drawdown  Data  Sources, 
Errors,  and  Analysis 


Prototype  data  on  navigation  effects  have  been  collected  on  the  Upper 
Mississippi,  Illinois,  and  Ohio  Rivers.  One  major  limitation  of  these  data  is 
that  only  a  few  points  were  collected  between  the  vessel  and  the  bank  line, 
primarily  in  the  near-bank  zone.  Consequently,  these  data  can  not  provide  a 
verification  of  the  velocity  distribution  shape  and  the  return  velocity 
magnitude  outside  the  near  bank  zone.  Physical  model  studies  of  reaches  on 
the  Upper  Mississippi  River  were  conducted  where  return  velocity  and 
drawdown  were  measured  from  shore  to  vessel  to  define  the  distribution  of 
return  velocity  and  drawdown.  The  physical  model  tests  simulated  a  reach  of 
the  Illinois  River  near  Kampsville  and  the  Mississippi  River  near  Clark’s 
Ferry  where  the  Illinois  State  Water  Survey  (ISWS)  (Bhowmik,  Soong,  and 
Xia  1993a,b)  conducted  prototype  measurements  of  tow-induced  return 
velocity  and  drawdown.  The  physical  model  was  verified  by  comparing 
model  and  prototype  results  for  six  tow  events  in  each  river  reach. 

Aereement  was  reached  between  model  and  prototype  return  velocity  and 
drawdown  by  reducing  the  physical  model  draft.  Details  of  the  Ksmpsville 
and  Clark’s  Ferry  physical  model  studies  and  the  data  used  herein  are  pro¬ 
vided  in  Maynord  and  Martin  (1996a,b). 

As  with  many  data  collection  efforts,  collection  of  return  velocity  and 
water  level  drawdown  data  is  not  an  easy  task  in  either  the  physical  model  or 
the  prototype.  Some  error  sources  are  specific  to  either  model  or  prototype  or 
applicable  to  both.  Errors  in  instrumentation  along  with  details  of  other 
potential  error  sources  in  the  physical  model  testing  are  discussed  in  the 
physical  model  reports. 

Instrument  error  in  the  ISWS  prototype  tests  must  be  considered  in  using 
prototype  data  because  of  the  environmental  conditions  found  in  the  prototype. 
Electromagnetic  velocity  meters  were  used  in  the  ISWS  prototype  studies  and 
can  be  subject  to  a  variety  of  problems.  Radio  interference  can  result  in 
highly  erratic  readings.  Drift  can  accumulate  around  the  sensor  head  and 
modify  readings  in  a  manner  that  is  not  obvious  to  the  person  conducting  the 
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tests.  Nevertheless,  the  electromagnetic  velocity  meter  is  the  best  available 
prototype  technique  for  this  type  of  measurement.  Wave  data  collected  in  the 
ISWS  tests  were  frequently  missing  due  to  instrument  malfunction. 

Tow  length  and  width  are  relatively  easy  to  determine  in  the  field  as  long 
as  the  tow  is  comprised  of  the  standard  10.7  by  59.5  (35-  by  l95-ft)  barges. 
When  chemical  barges  are  present,  their  nonuniformity  of  length  and  width 
can  be  a  problem.  Tow  length  is  important  because  the  length  and  time  of 
passage  past  a  given  point  are  often  used  to  determine  the  tow  speed  in  proto¬ 
type  tests.  Tow  draft  may  also  contribute  to  the  data  scatter.  Tows  are 
assumed  to  draft  2.7  m  (9  ft)  if  loaded  and  0.6  m  (2  ft)  if  unloaded  unless 
contact  with  the  towboat  captain  indicates  otherwise.  Mixed  tows  further 
complicate  the  issue. 

Also  the  vertical  location  of  velocity  measurements  is  considered.  Varia¬ 
tions  caused  by  this  effect  are  not  as  significant  as  it  might  first  appear  since 
tow-induced  return  velocity  profiles  are  generally  uniform.  The  boundary 
layer  resulting  from  tow  passage  does  not  have  sufficient  time  to  completely 
develop  into  the  typical  velocity  profile  found  in  a  fully  developed  open  chan¬ 
nel  flow.  Therefore,  velocity  changes  near  the  surface  should  be  uniform 
down  to  the  point  where  the  boundary  layer  has  developed.  Results  from  the 
ISWS  tests  provided  in  Maynord  and  Martin  (1996a)  suggest  that  only  meters 
1(X)1  and  998  (Figure  2)  in  the  trip  1  data  should  not  be  used  in  the  compari¬ 
son  with  the  1991  analytical  method  because  of  the  boundary  layer  influence. 

One  of  the  greatest  causes  of  data  variability  is  determining  the  tow  impact 
from  the  time-history  of  the  parameter  of  interest.  The  difficulty  of  this  task 
increases  as  the  magnitude  of  the  tow  impact  decreases,  because  it  becomes 
more  difficuit  to  distinguish  the  tow  impact  from  the  natural  variations  present 
in  the  river.  The  first  task  in  defining  the  tow  impact  is  defining  the  ambient 
velocity  prior  to  tow  passage  which  can  be  a  problem  as  there  are  some  long 
period  variations  in  the  ambient  velocity  due  to  several  factors  such  as 
changes  in  the  regulating  gates  at  the  locks  and  dams.  The  next  task  is  defin¬ 
ing  the  maximum  deviation  from  ambient  conditions.  In  the  ISWS  tests,  the 
maximum  (or  minimum)  value  from  an  11-sec  moving  average  was  selected 
during  the  time  the  tow  was  adjacent  to  the  meter.  This  smoothing  technique 
eliminated  the  normal  variations  present  in  turbulent  flow  while  also  defining 
the  maximum  value.  The  difference  between  the  maximum  value  and  the 
ambient  was  the  tow  impact. 


Chapter  3  Return  Velocity  and  Drawdown  Data  Sources,  Errors,  and  Analysis 


Depth  <n> 


Notei  Meters  998,  999,  and  1000  are  nounted  on  ^  platforn 

ot  0.31,  1.22,  and  2.44  n  above  the  botton,  respectively. 
The  velocity  neters  are  two  dinensional  elec tronagne tic 
neters.  Details  of  the  neter  descriptions  can  be  found 
in  Bhownik,  Soong,  and  Xla  <1993a>. 


Figure  2.  Cross  section  of  the  Illinois  River  at  the  Kampsville  site  for  trip  1 
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4  1991  and  1995  Analytical 

Methods 


Return  Current  Comparisons  with  1991  Analytical 
Method 

Tow  speed  used  in  application  of  the  analytical  methods  was  equal  to  the 
sum  of  vessel  speed  relative  to  the  ground  and  the  ambient  velocity.  Ambient 
velocity  was  positive  for  unbound  tows  and  negative  for  downbound  tows. 
Data  were  limited  to  tows  with  a  length  at  least  40  percent  of  the  channel 
width  to  eliminate  the  tow  length  effects  measured  in  physical  model  tests  of 
the  Clark’s  Ferry  reach  being  studied  on  the  Mississippi  River  (Maynord  and 
Martin  1996b). 

A  scatterplot  of  observed  return  velocity'  versus  the  predicted  return 
velocity  using  the  1991  analytical  method  (Figure  3)  was  developed  from  the 
Kampsville  physical  model  data.  The  1991  analytical  method  under  predicts 
return  velocity  observed  in  the  physical  model  for  most  of  the  tests. 

Comparison  of  the  1991  analytical  method  to  the  ISWS  data  required  deter¬ 
mining  the  draft,  beam,  and  length  to  use  for  some  unusual  tow  conflgurations 
having  a  mixture  of  loaded,  unloaded,  and  partially  loaded  barges.  The 
approach  for  this  report  selected  the  beam,  draft,  and  length  of  the  section 
having  the  maximum  cross-sectional  area. 

A  scatterplot  (Figure  4)  of  observed  prototype  return  velocity  versus  the 
predicted  return  velocity  from  the  1991  analytical  method  was  developed  from 
ISWS  Kampsville  prototype  test  data.  As  in  the  physical  model,  prototype 
data  were  limited  to  tows  with  a  length  greater  than  40  percent  of  the  chwnel 
width.  N  less  than  or  equal  to  51  was  used  to  agree  with  the  limits  of  the 
physical  model  data.  Meters  998  and  1001  from  trip  1  were  excluded  because 
their  close  proximity  to  the  boundary  (0.31  m)  led  to  boundary  layer  effects. 
As  in  the  physical  model,  the  1991  analytical  method  under  predicts  observed 
prototype  return  velocity  for  most  tests. 
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OBSERVED  RETURN  VELOCITY.  M/SEC 


Figure  3,  Observed  physical  model  return  velocity  versus  computed  return 
velocity  using  1991  analytical  method 


1995  Analytical  Method 

Based  on  the  underprediction  of  both  physical  model  and  prototype  results, 
the  1991  ana'j'tical  method  must  be  reexamined.  The  physical  model  data 
from  Kampsville  and  Clark’s  Ferry  provided  an  opportunity  to  evaluate  each 
pari  of  the  1991  method.  The  maximum  return  velocity  'V,(Y)  during  the  tow 
event  for  each  meter  is  the  value  extracted  from  the  model  and  prototype  data. 
The  V,(Y)  data  from  the  eight  physical  model  meters  from  the  Kampsville 
and  Clark’s  Ferry  physical  model  data  were  averaged  to  determine  the  return 
velocity  representative  of  the  entire  cross  section  V,  Since  V,(Y)  is  the 
maximum  return  velocity  during  the  tow  event  at  Y,  V,  is  a  cros.s-sectional 
average  of  these  maximums.  The  ratio  Vr/Schijf  Vr  (from  Eq’s  1  and  2)  was 
compared  for  upbound  versus  downbound  vessels  to  determine  if  the  average 
channel  velocity  is  the  correct  velocity  to  use  in  determining  the  velocity  of 
the  tow  relative  to  the  water.  A  likely  alternative  is  the  velocity  where  the 
tow  is  operating  which  is  generally  greater  than  the  average  channel  velocity. 
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Figure  4.  Observed  prototype  return  velocity  versus  computed  return 
velocity  using  1991  analytical  method 

The  average  Vr/Schijf  Vr  for  upbound  tows  was  higher  than  the  average 
ratio  for  downbound  tows  when  using  the  average  channel  velocity  to  deter¬ 
mine  the  vessel  speed  through  the  water.  When  the  average  channel  velocity 
was  increased  by  20  percent  and  then  used  to  determine  the  speed  through 
the  water,  upbound  and  downbound  tows  gave  the  same  average  ratio.  This 
finding  tends  to  confirm  the  use  of  the  velocity  where  the  tow  operates  when 
determining  the  tow  speed  through  the  water.  The  20-percent  increase  is 
consistent  with  the  actual  velocity  that  occurs  in  the  deeper  areas  where  the 
tow  is  operating 

The  Schijf  return  velocity,  calculated  using  Equations  I  and  2  with  the 
20-percent  increase  in  Vavg,  under  predicted  the  observed  cross-sectional 
average  of  the  maximum  Vr  from  the  physical  model.  To  correct  this  under¬ 
prediction,  an  effective  draft  and  beam  were  determined  by  adding  the  dis¬ 
placement  thickness  to  the  draft  and  twice  the  displacement  thickness  to  the 
beam  of  the  vessel  as  recommended  by  Gates  and  Herbich  (1977).  Displace¬ 
ment  thickness  is  determined  from  the  PrandtI-Schlicting  skin  friction 
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equation  for  a  smooth  flat  plate  at  zero  incidence  (Schlicting  1968).  Dis< 
placement  thickness  6,  is  determined  from  the  skin  fnction  equation 
according  to 


.  0.292  L 


(12) 


where 


L  -  total  barge  length 
Rl  =  V,,,  Uv 
^^,,  =  V,+V^,+SchijfVr 
Vg  -  vessel  speed  over  the  ground 

^tmb  ~  positive  for  upbound  tows  and  negative  for  downbound  tows 
V  =  kinematic  viscosity  c  f  water 

While  this  significantly  improved  the  comparison  between  observed  average 
K,  and  P',(SchijO.  the  error  between  the  two  was  a  function  of  the  ratio  of 
vessel  speed  to  the  critical  or  limiting  speed  (Vi)  from  Equations  1  and  2. 
Dependence  with  is  consistent  with  the  findings  of  Schijf  (1949).  The 
observed  physical  model  data  were  used  to  develop  a  correction  factor  based 
on  V/Vi  which  provided  good  agreement  between  observed  average  F,  and 
Schijf  Vf.  Schijf  V,  is  calculated  using  the  effective  draft  and  beam  in  the 
Schijf  equation.  Average  V,  is  calculated  from 


V,  =  V^iSchijf)  (1.9-1.29V/V^) 


(13) 


The  computed  F,  from  Equation  13  is  not  allowed  to  be  less  than  the  return 
velocity  from  the  Schijf  equation  and  V/VL  should  be  from  0.35  to  0.9.  The 
comparison  between  observed  F,  from  the  Kampsville  and  Clark’s  Ferry 
physical  models  and  computed  V,  using  Equation  13  is  shown  in  Figure  S. 

S.'  d,  is  computed  using  F,  from  equations  1  and  2,  the  solution  does  not 
req..ire  iteration  with  Equation  13. 

Hie  next  step  is  evaluating  the  relationship  in  Equation  8  that  proportions 
the  return  velocity  on  each  side  of  the  vessel.  The  return  velocity  for  each 
side  of  the  vessel  F„  was  determined  for  each  test  using  the  Kampsville  and 
Clark’s  Ferry  physical  model  data.  F„  was  normalized  by  the  observed 
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Figure  5.  Observed  physical  model  average  return  velocity  versus  computed 
Schijf  return  velocity 


average  V,  (which  represents  the  entire  cross  section)  and  tested  against 
various  parameters  such  as  SKEW  used  in  Equation  8.  The  best  fit  was 
found  using  the  relation  between  distance  from  vessel  to  shoreline  and 
total  channel  width  For  B,tj/B,^^  from  0  to  0.5 


=  1.65-1.3^2* 

V. 


For  from  0.5  to  1.0 
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(15) 


X 


1.33-0.7^ 


Equations  14  and  15  are  plotted  against  observed  data  from  the  Kampsville 
and  Clark’s  Ferry  physical  models  in  Figure  6  along  with  rectangular  flume 
data  from  a  previously  unpublished  data  set. 


Figure  6.  versus  using  physical  model  data 


The  third  step  is  defining  a'  “  V„JV„  which  defines  the  uniformity  of  the 
return  velocity  distribution  between  vessel  and  shoreline.  Previous  attempts 
at  this  relation  using  physical  model  data  were  only  approximate  because 
scale  effects  were  not  eliminated  from  the  data  as  has  been  done  in  the 
Maynord  and  Martin  (i996a,b)  physical  model  data.  Only  the  physical 
model  data  from  the  Kampsville  model  were  used  in  the  analysis  of  a' 
because  of  the  likely  influence  of  dikes  from  the  Clark's  Ferry  model.  The 
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physical  model  data  were  used  to  determine  V„„  at  (Y-B)/(B,jj,-B)  =  0. 

Various  parameters  were  tested  against  a',  and  provided  the  best 
agreement  with  observed  data.  V„JV„  is  plotted  against  in  Figure  7  and 

is  described  by  the  equation  # 


a 


=  0J5NJ'* 


(16) 
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Figure  7.  versus  using  physical  model  data 


The  final  step  was  determining  C  in  Equation  10  which  describes  the 
shape  of  the  return  velocity  distribution.  One  requirement  of  C  is  that  it 
must  equal  zero  for  a'  =  1 ,  The  C  that  best  defines  the  shape  of  the  return 
velocity  distribution  in  the  physical  model  date  from  Kampsville  and  Clark’s 
Ferry  is  defined  by 
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The  revised  equations  presented  above  will  be  referred  to  hereafter  as  the 
1995  analytical  method.  A  scaherplot  of  the  Kampsville  physical  model  data 
versus  the  1995  analytical  method  (Figure  8)  shows  an  improved  comparison 
of  observed  versus  predicted  return  velocity.  Individual  test  plots  for  the 
Kampsville  physical  model  are  shown  in  Figures  9  to  13.  The  six  plots  on 
each  figure  for  individual  tests  having  varying  speed,  position,  upbound 
versus  downbound,  etc.  In  the  test  number,  the  K  stands  for  Kampsville  and 
L  or  H  for  low  or  high  flow.  R  or  L  stand  for  right  or  left  of  the  channel 
thalweg  and  U  or  D  stand  for  upbound  or  downbound.  If  neither  R  or  L 
appear  before  U  or  D,  the  tow  is  on  the  channel  thalweg.  The  last  three 
numbers  give  the  vessel  speed  relative  to  the  ground  in  model  meters  per  . 
second  which  can  be  converted  to  prototype  values  for  the  1  ;25  Kampsville 
model  by  multiplying  by  25'"=  5.  Details  of  each  test  are  given  in  Table  1. 
On  some  of  the  physical  model  tests,  the  exponential  type  equation  given  by 
Equation  10  underpredicts  the  return  velocity  at  the  bank.  Other  data  sets 
will  be  examined  to  see  if  this  trend  is  repeated.  A  scatterplot  of  the  Clark’s 
Ferry  physical  model  data  versus  the  1995  analytical  model  is  shown  in  Fig¬ 
ure  14  for  N  less  than  52.  Figure  15  provides  a  scatterplot  of  tests  con¬ 
ducted  with  pool  572.7  in  the  Clark’s  Ferry  section  which  resulted  in  N  of 
about  85.  The  tests  in  Figure  15  addressed  the  open  river  portion  of  the 
Upper  Mississippi  River  where  large  channel  cross  sections  result  in  large 
blockage  ratios.  The  1995  analytical  method,  derived  from  physical  model 
data,  was  then  evaluated  using  the  ISWS  data.  Scatterplots  of  prototype  data 
from  Kampsville  and  Clark’s  Ferry  are  shown  in  Figures  16  and  17,  respec¬ 
tively.  Individual  plots  for  the  ISWS  prototype  data  at  Kampsville  are 
shown  in  Figures  18  to  21.  Each  of  the  24  plots  in  Figures  18  to  21  repre¬ 
sent  an  individual  prototype  tow  with  the  towboat  name  given  in  the  upper 
right-hand  corner.  Plots  where  there  are  two  or  more  points  at  the  same 
distance  from  the  tow  are  locations  where  velocity  meters  were  placed  at 
different  vertical  positions.  Details  of  the  prototype  tests  are  provided  in  the 
Table  2  report.  Note  that  the  Kampsville  prototype  data  do  not  show 
evidence  of  an  increase  near  the  bank. 
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KAMPSVILLE  PHYSICAL  MODEL  POOL  418.0 
aODOD  KAMPSVILLE  physical  model  pool  419.4 
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Figure  8.  Observed  Kampsville  physical  model  return  velocity  versus  computed  return 
velocity  using  1995  analytical  method 
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Figure  9.  Observed  Kampsville  physical  model  return  velocity  versus  computed  return  velocity 
using  1995  analytical  method  for  individual  tests:  KLU335,  488,  640,  354  and 
KLD506  and  659  (Continued) 
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VELOCITY 


e.  Test  KLD506  F.  Test  KLD659 

Figure  9.  (Concluded) 
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a.  Test  KLRU49 


b.  Test  KLRD49 


c.  Test  KLLU49  d.  Test  KLLD51 

Figure  1 0.  Observed  Kampsville  physical  model  return  velocity  versus  computed  return 
velocity  using  1 995  analytical  method  for  individual  tests;  KLRU49,  KLRD49, 
KLLU49,  KLLD51,  KL1U46,  and  KL1U61  (Continued) 
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Figure  10.  (Concluded) 
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a.  TestKL1U76  b.  Test  KLEU49 


c.  Test  KLEU67  d.  Test  KHEU38 

Fig- '  b  • .  Observed  Kampsville  physical  model  return  velocity  versus  computed  return 
velocity  using  1995  analytical  method  for  individual  tests:  KL1U76,  KLEU49, 
KLEU67,  KHEU38,  KHEU56,  and  KHOU38  (Continued) 
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VELOCITY 
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Figure  1 1 .  (Concluded) 
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a,  Test  KHOU53 


b.  Test  KHOD50 


c.  Test  KH0D66  d.  Test  KHRU38 

Figure  1 2.  Observed  Kampsville  physical  model  return  velocity  versus  computed  return 

velocity  using  1995  analytical  method  for  individual  tests:  KHOU53,  KHOD50, 
KHOD66,  KHRU38,  KHRD66,  and  KHLU38  (Continued) 
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e.  Test  KHRD66  f.  Test  KHLU38 

Figure  12.  (Concluded) 
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Figure  13.  Observed  Kampsville  physical  model  return  velocity  versus  computed  return 

velocity  using  1995  analytical  method  for  individual  tests:  KHLD66,  KH0U27, 
KH0064,  LU38C2,  and  LD58Q2  (Continued) 
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PREDICTED  RETURN  VELOCITY.  M/SEC 


1.0 


*****  CLARK'S  FERRY  PHYSICAL  MODEL  POOL  546.0 
DDDOO CLARK'S  FERRY  PHYSICAL  MODEL  POOL  651.5 


OBSERVED  RETURN  VELOCITY.  M/SEC 


Figure  14.  Observed  Clark's  Ferry  physical  model  return  velocity  versus  computed  return 
velocity  using  1995  analytical  method  for  N  <  52 
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PREDICTED  RETURN  VELOCITY.  M/SEC 


Figure  15.  Observed  Clark's  Ferry  physical  model  return  velocity  versus  computed  return 
velocity  using  1 995  analytical  method  for  N  =  85 
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*♦**♦  KAMFSVILLE  TRIP  1  PROTOTYPE  DATA 
□□□□□  KAMPSVILLE  TRIP  2  PROTOTYPE  DATA 


OBSERVED  RETURN  VELOCITY,  M/SEC 


Figure  1 6.  Observed  Kampsville  prototype  return  velocity  versus  computed  return 
velocity  using  1 995  analyticai  method 
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PREDICTED  RETURN  VELOCITY,  M/SEC 


*****  CLARK'S  FERRY  TRIP  1  PROTOTYPE  DATA 


OBSERVED  RETURN  VELOCITY,  M/SEC 


Figure  17.  Observed  Clerk's  Ferry  prototype  return  velocity  versus  computed  return 
velocity  using  1 995  analytical  method 
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VQ-OCl 


a.  MA/  Mr.  Aldo  b.  MA/  Floyd  Btaske 
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c.  MA^  Sugarl  (13)  d.  MA/  W.  C.  Norman  (13) 

Figure  1 8.  Observed  Kampsville  prototype  return  velocity  versus  computed  return  velocity 
using  1995  analytical  method  for  individual  tests:  M/V  Mr.  Abdo,  Floyd  Blaske, 
Sugar  (13),  W.  C.  Norman  (13),  ContiKarla,  and  Rambler  (Continued) 
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VELOCITY 


a.  M/VMIawrce(14)  b.  MA^  ChaLehman 


c.  Jeffboat 


d.  M/V  Ardyce  Rondl 


Figure  1 9.  Observed  Kampsville  prototype  return  velocity  versus  computed  return  velocity 
using  1 995  analytical  method  for  individual  tests:  MA/  Mlawrce,  ChaLehman, 
Jeffboat,  Ardv-e  Randl,  Mr.  Paul  (15),  and  Marget  0  (15)  (Continued) 
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Figure  19.  (Concluded) 


(0 


Chapter  4  1991  and  1995  Analytical 


M/V  Dixie  Patrit 


M/V  Orleanian 


FiQure  20.  Observed  Kampsville  proiotype  return  velocity  versus  computed  return  velocity 

using  1995  analytical  method  for  Individual  tests:  MA/  Mr.  Lawre  (15),  Al  Smith, 
Dixie  Patrit,  Orleanian,  Pat  Breen,  and  Dixie  Expre  1  (Continued) 
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Figure  20.  (Concluded) 
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•  •••*  ICArt»SVIU-E  PROTOTYPE  DATA 
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Figure  21.  Observed  Kampsvilie  prototype  return  velocity  versus  computed  return  velocity 

using  1995  analytical  method  for  individual  tests:  MA/  Night  1401,  Irving  Crown, 
Thurston  Men,  Olmstead,  Jack  D.  Wofford,  and  Dixie  Expre  2  (Continued) 
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e.  M/V  J&Lk  D.  Woffod  f.  M/V  Dixie  Expre  2 

Figure  21.  {Concluded) 
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Table  2  (Concluded) 
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i-ri' . 

1 

■ 

'• 

1  ■  , 


! 


1 
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5 


Development  of  1 995 
Analytical  Method  for 
Prediction  of  Water-Level 
Drawdown 


To  compute  water-level  drawdown  at  position  Y  from  the  center  line  of  the 
vessel,  die  average  return  velocity,  as  d«^termined  from  Equation  13,  was  used 
in  Equation  2  to  determine  the  average  o'rawdown  z  for  the  section.  The 
difficulty  at  this  point  is  that  the  numbtr  af  drawdown  distributions  in  the 
physical  model  was  limited.  Therefore,  return  velocity  distributions  were 
us^  to  determine  how  the  drawdown  varied  on  either  side  of  the  tow  and 
across  the  section.  The  drawdown  on  either  side  of  the  tow  was  determined 
from  the  same  equations  used  for  return  velocity  in  Equations  14  and  IS. 

For  from  0  to  0.5 

h  =  1.65-1.3^  (18) 

2  B^. 


For  B^B^  from  0.5  to  1.0 

h.  =  1.35-OjE^  (19) 

Z 


The  maximum  drawdown  at  (Y-B)I{B^,-B)  =®  0  is 


Ch«pt«r  G  Davaiopfnent  of  1 995  Analytical  Method  for  Prediction 
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(20) 


where  is  a  function  of  a'  defined  by  Equation  16. 
The  distribution  of  Z(Y)  is  defined  as 


^  =  «p  [C(-^)l 


where  Z(Y)  is  the  maximum  water-level  drawdown  during  the  tow  event  at 
distance  Y  from  the  vessel  center  line  and  C  is  defined  by  Equation  17. 
is  the  primary  factor  that  establishes  the  slope  of  the  drawdown  distribution. 
Based  on  the  Kampsville  and  Clark’s  Ferry  physical  model  drawdown  data 
and  detailed  drawdown  measurements  for  the  WC  Norman  tow  in  the 
Kampsville  physical  model,  is  defined  by 


Scatterplots  of  observed  physical  model  drawdown  from  Kampsville  and 
Clark’s  Ferry  versus  computed  drawdown  using  the  1995  analytical  drawdown 
method  are  shown  in  Figures  22  and  23,  respectively.  Comparison  of  the 
analytical  method  with  detailed  drawdown  data  from  the  WC  Norman  tow  in 
the  Kampsville  physical  model  is  shown  in  Table  3.  The  physical  model- 
based  1995  Analytical  Drawdown  method  was  then  tested  for  agreement  with 
the  prototype  data  which  was  taken  at  a  single-wave  gauge  at  the  Kampsville 
site  on  the  Illinois  River.  The  wave  data  from  ISWS  were  collected  at  10 
samples/sec.  To  be  consistent  with  the  return  velocity  analysis,  an  11-sec 
moving  average  was  used  to  smooth  the  wave  data  before  extracting  the 
ambient  and  changed  water  level  from  the  record.  Drawdown  measurements 
based  on  only  a  staff  gage  reading  were  not  used  in  this  comparison.  A 
scatterplot  of  drawdown  for  the  ISWS  prototype  data  is  shown  in  Figure  24. 
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PREDICTED  DRAUJOOUN, 


*****  KAUPSVILLE  physical  model  pool  418.0 
□ODOD  KAMPSVILLE  PHYSICAL  MODEL  POOL  419,4 
AAA4A  KAMPSVILLE  PHYSICAL  MODEL  POOL  427.0 


Figure  22.  Observed  Kampsville  physical  model  drawdown  versus  computed  drawdown 
using  1995  analytical  method 
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CLARK'S  FERRY  PHYSICAL  MODEL  POOL  546.0 
□□□DO CLARK'S  FERRY  PHYSICAL  MODEL  POOL  551.5 
AAZkA^aARK'S  FERRY  PHYSICAL  MODEL  POOL  572.7 


OBSERVED  DRAWDOWN,  M 


Figure  23.  Observed  Clark's  Ferry  physical  model  drawdown  versus  computed  drawdown 
using  1995  analytical  method 
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Chapter  5  Development  of  199b  Aniaytical  Method  for  Prediction 


1  Table  3 

Observed  Versus  Computed  Drawdown  for  WC  Norman  Test 


Dittance  from  tow,  m 

-162 

-100 

-41 

34 

70 

106 

124 

Obaarvad  phya  model,  cm' 

S.4 

6.5 

6.8 

7,7 

6.8 

6.9 

- 

Obaerved  Prototype,  cm' 

- 

- 

- 

- 

- 

- 

7.2 

Computed  drawdown,  cm' 

6.4 

6.1 

6.S 

B.4 

7.7 

7.2 

6.9 

Note:  (-1  No  vohit  wat  avaHiUa. 
'  Thaaa  vakiaa  art  drawdown. 


*****  kampsville  prototype  trip  1 


OBSERVED  DRAUIDOWN,  fl 


Figure  24.  Observed  Kampsville  prototype  drawdown  versus  computed  drawdown  using 
1995  analytical  method 
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6  Comparison  of  1995 
Analytical  Method  with 
Independent  Field  Data 


Since  the  Kaiq>sville  and  Clark’s  Ferry  field  data  were  used  to  verify  the 
Kan^ttville  and  Clark’s  Ferry  physical  m^els  which  were  then  used  to  pro¬ 
duce  data  used  in  the  development  of  the  1995  method,  comparison  of  the 
analytical  method  to  the  same  field  data  is  not  independent.  Two  field  data 
sets  were  used  to  provide  an  independent  conqiarison  with  the  1995  analytical 
method.  The  first  data  set  was  return  velocity  data  taken  by  ESE  (1981)  at 
one  section  on  the  Mississippi  River  and  one  on  the  Illinois  River.  A 
scatteiplot  of  all  ESE  return  velocity  data  meeting  the  limitations  of  the  1995 
method  is  shown  in  Figure  25.  The  second  data  set  comes  from  data  collected 
at  four  sites  on  the  Ohio  River  by  the  U.S.  Army  Engineer  District, 

Louisville,  whose  cross  sections  are  shown  in  Figures  26  through  29  and  data 
on  Tables  4  through  7.  Scatterplots  for  all  four  sections  are  shown  in 
Figures  30  and  31  for  return  velocity  and  drawdown,  respectively.  Scatter  is 
large,  but  the  trend  around  the  line  of  perfect  agreement  is  correct  for  all 
three  scatterplots. 
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Chaptar  6  Compartaon  of  1 995  Analytical  Method  with  Independent  Field  Data 


**»«c*ESE  MISSISSIPPI  RIVER 


OBSERVED  RETURN  VELOCITY.  M/SEC 


Figure  25.  Observed  ESE  prototype  return  velocity  versus  computed  return  velocity  using 
1995  analytical  method 
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LErr  BAM< 


Figure  26.  Cross  section  for  ORL  tests  89192-89196,  French  Island  (Factor 
for  converting  feet  to  metric  is  3.048  and  square  feet  to 
square  meters  is  0.0929) 
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square  meters  is  0.0929) 
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Figure  28.  Cross  section  for  ORL  tests  90172*90179  (Factor  for  converting 
feet  to  meters  is  3.048  and  square  feet  to  square  meters  is 
0.0929) 
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Figure  29.  Cross  section  for  ORL  tests  90191-90195  (Factor  for  converting 
feet  to  meters  is  3.048  and  square  feet  to  square  meters  is 
0.0929) 
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Tabto4 

ORL  Twt  89192-89196 

89192003  108 


9999  mains  data  not  racordad  or  matar  malfunction. 
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Table  4  (Concluded) 


1390 

0.65 

0 

10 

6 

Bi 

1445 

2.6 

16 

iS 

89195008 

105 

1170 

9 

999 

8.8 

1.3 

1.0 

14100 

465 

36360 

1465 

1095 

0.5 

0 

18 

4 

1390 

0 

10 

6 

m 

1445 

0.75 

16 

m 

B91 96001 

175 

i120 

2 

999 

10.5 

1.8 

-1.0 

29700 

965 

36360 

1466 

1095 

0.8 

0 

18 

4 

1390 

0.55 

0 

10 

6 

1445 

2 

16 

89196002 

108 

1082 

2 

999 

16.7 

1.8 

-1.0 

19200 

615 

36360 

1465 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

1445 

0 

16 

89196003 

54 

720 

9 

999 

13.5 

1.8 

1.0 

25800 

615 

36360 

1465 

1095 

mm 

0 

18 

4 

1390 

0.3 

0 

10 

6 

1445 

2 

16 

89196096 

105 

1116 

9 

999 

■■ 

12.5 

1.8 

•1.0 

15000 

490 

36360 

1465 

1390 

0.5 

im 

10 

6 

9993 

9999 

9999 

9999 

9999 

1445 

1.5 

16 

89196007 

64 

1129 

9 

999 

13.9 

1.8 

1.0 

15400 

516 

36360 

1465 

1095 

0.45 

0 

18 

4 

m 

1390 

0.55 

0 

10 

6 

m 

1445 

2.25 

16 

Note:  To  convert  to  metric,  multiply  by  these  factors: 
feet  X  0.3048  to  obtain  meters 
feet'  X  0.0929  to  obtain  square  meters 
horsepower  x  9809.5  to  obtain  watts 
inches  x  0.0254  to  obtain  meters 
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Chapter  6  Comparison  of  1 995  Analytical  Method  with  Independent  Field  Data 


Table  5 

ORL  Test  89198-89203 


89198002 

54 

1000 

!■ 

999 

11 

1.1 

1.0 

12200 

495 

20460 

1020  11 

76 

0.25 

0 

13 

8 

20 

1.5 

4 

• 

B 

89 198003 

105 

1139 

9 

999 

12.6 

1.1 

-1.0 

13500 

545 

20460 

1020 

76 

0.65 

0 

13 

8 

20 

4.5 

4 

89198004 

175 

1113 

2 

999 

11.4 

n 

-1.0 

12900 

520 

20460 

1020 

76 

mm 

0 

13 

8 

20 

2 

4 

89198005 

194 

1115 

9 

999 

7.25 

1.1 

-1,0 

12900 

520 

20460 

1020 

76 

0.5 

0 

13 

8 

,  . 

20 

1.5 

4 

89199001 

105 

1125 

9 

999 

12.9 

-1.0 

13600 

365 

20460 

1020 

76 

0.6 

0 

13 

8 

20 

5 

4 

89199002 

35 

269 

2 

999 

15.5 

-1.0 

470 

20460 

1020 

9999 

9999 

9999 

9999 

9999 

20 

0 

m 

89199003 

245 

1120 

9 

999 

7.52 

IH 

1.0 

10300 

430 

20460 

1020 

76 

1.06 

tm 

13 

8 

20 

4 

4 

89199004 

105 

1065 

8 

999 

9.05 

ilM 

-1.0 

11700 

480 

20460 

1020 

76 

EBi 

0 

13 

8 

20 

1.5 

4 

89200002 

70 

9 

999 

9.5 

1 

-1.0 

15800 

650 

20460 

1020 

76 

0 

13 

8 

20 

1.5 

4 

69200003 

105 

1155 

9 

999 

12.2 

1 

1.0 

14100 

570 

20460  ' 

1020 

76 

0.55 

0 

13 

8 

20 

2.25 

4 

89201001 

105 

1139 

2 

999 

11 

1.8 

1.8 

14500 

595 

20460 

1020 

9999  means  data  not  recorded  or  meter  malfunction. 
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Chapter  6  Conip8''isc'n  of  1995  Analytical  Method  with  Independent  Reid  Data 


Table  5  (Concluded) 


89202008 

106 

892 

9 

999 

6.21 

1.8 

1.0 

12900 

520 

20460 

1020 

9999 

9999 

9999 

9999 

9999 

20 

1 

4 

89203001 

34 

148 

9 

999 

12.8 

1.8 

13100 

530 

20460 

1020 

9999 

9999 

9999 

9999 

9999 

n 

2C 

0 

4 

89203002 

105 

1142 

9 

999 

10.2 

1.8 

-1,0 

11700 

480 

20460 

1020 

9999 

9999 

9999 

9999 

9999 

20 

2 

4 

89203003 

70 

892 

am 

999 

9.1 

1.8 

1,0 

13400 

540 

20460 

1020 

9999 

9999 

9999 

9999 

9999 

20 

1.5 

4 

89203004 

105 

1070 

9 

999 

7.04 

1.8 

1.0 

12700 

480 

20460 

1020 

9999 

9999 

9999 

9999 

9999 

20 

1.5 

wm 

89203005 

105 

1113 

9 

999 

9.8 

1.8 

-1,0 

10200 

420 

20460 

1020 

9999 

9999 

9999 

9999 

9999 

20 

1 

BH 

Note:  To  convert  to  metric,  multiply  by  these  factors: 
feet  X  0.3048  to  obtain  meters 
feet^  X  0.0929  to  obtain  square  meters 
horsepower  x  9809.5  to  obtain  watts 
inches  x  0.0254  to  obtain  meters 
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Tabl«  6 

ORL  Test  90172-90179 

901 72004 

105 

687 

9 

999 

7.61 

0.9 

1.0 

7500 

350 

28626 

1100 

130 

0.25 

0 

23 

6 

960 

0.16 

0 

25 

6 

9999 

9999 

9999 

90172006 

105 

705 

9 

999 

9.13 

0.9 

1.0 

4200 

225 

28626 

1100 

960 

0.25 

0 

25 

6 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

90172009 

108 

870 

9 

999 

12.3 

0.9 

1.0 

1000 

100 

28626 

1100 

130 

0.75 

0 

23 

6 

960 

0.4 

0 

25 

6 

20 

S 

4 

90173002 

105 

1113 

2 

999 

8.55 

0.7 

1.0 

7500 

350 

28626 

1100 

130 

0.15 

0 

23 

6 

960 

0.1 

0 

25 

6 

20 

0.5 

4 

901 7400 J 

105 

1113 

9 

999 

11.7 

0.9 

1.0 

7500 

350 

28626 

1100 

130 

0.7 

0 

23 

6 

960 

0.4 

0 

25 

6 

20 

4.5 

■ 

4 

901 74004 

106 

1139 

9 

999 

12.1 

0.9 

•1.0 

6700 

325 

28626 

1100 

130 

0.4 

0 

23 

6 

960 

0.35 

0 

25 

6 

20 

3.5 

4 

90174005 

105 

1080 

9 

999 

9,97 

0.9 

-1.0 

8900 

400 

28626 

1100 

130 

0.4 

0 

23 

6 

960 

0.4 

0 

25 

6 

20 

2 

4 

90175002 

105 

1115 

9 

999 

7.9 

1 

1.0 

8900 

400 

28626 

1100 

130 

0.25 

0 

23 

6 

9999  means  data  not 

recorded 

or  meter 

malfunct 

on. 

(Continued) 
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Table  6  (Concluded) 


960 

0.2 

0 

23 

6 

20 

2 

4 

90176001 

105 

1125 

9 

999 

9.92 

0.7 

1.0 

7500 

350 

28626 

1100 

130 

0.4 

0 

23 

6 

960 

0.3 

0 

25 

6 

20 

1.5 

4 

90176002 

105 

1105 

9 

999 

11.2 

0.7 

-1.0 

8800 

400 

28626 

1100 

130 

0.4 

0 

23 

6 

960 

0.3 

0 

25 

6 

20 

1.6 

4 

90176003 

35 

687 

9 

999 

10.8 

0.7 

1.0 

4800 

250 

28626 

1100 

130 

0.15 

0 

23 

6 

960 

0 

0 

25 

6 

20 

2 

4 

901 77002 

105 

1115 

9 

999 

11.5 

0.7 

-1.0 

11600 

500 

28626 

1100 

130 

0.35 

0 

23 

6 

960 

0.25 

0 

25 

6 

20 

2.5 

4 

90179005 

108 

1085 

9 

999 

12.1 

0.6 

-1.0 

6800 

330 

28626 

1100 

960 

0.3 

0 

25 

6 

130 

0.4 

0 

23 

6 

20 

3 

4 

90179006 

70 

738 

9 

999 

18.6 

0.6 

1.0 

7500 

350 

28626 

1100 

130 

0.4 

0 

23 

6 

960 

0.15 

0 

25 

6 

20 

3 

4 

Note:  To  convert  to 
feet  X  0.3048 
feet’  X  0.0929 
horsepower  x 
inches  x  0.025 

netric,  m 
to  obtain 
to  obtair 
3809.5  t( 
4  to  obtE 

jitiply  by  these  factors: 
meters 

square  meters 
obtain  watts 
in  meters 
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Chapter  6  Comparison  of  1995  Analytical  Method  with  Independent  Field  Data 


1  Table  7  (Continued) 

. .  ~l 

901 92003 

105 

1139 

2 

999 

8.9 

1.1 

1.0 

14400 

1100 

26346 

1850 

350 

0.2 

0 

13 

4 

125 

0.3 

0 

12 

4 

1530 

0.45 

0 

18 

4 

1725 

0.3 

0 

17 

4 

1830 

2.5 

10 

90192004 

105 

760 

9 

999 

13.2 

1.1 

•1.0 

14400 

1100 

26346 

1850 

350 

0.1 

0 

13 

4 

125 

0.25 

0 

12 

4 

1530 

0,75 

0 

18 

4 

1725 

0.4 

0 

17 

4 

1830 

3.25 

10 

90193001 

SO 

914 

9 

999 

10.4 

1.1 

1.0 

13375 

975 

28196 

1850 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

350 

0.15 

0 

13 

4 

1 

1830 

1.5 

10 

90193002 

106 

700 

2 

999 

13.1 

1.1 

-1.0 

13800 

1000 

28196 

1850 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

350 

0 

0 

13 

4 

125 

0 

0 

12 

4 

1830 

0.75 

10 

90193003 

105 

1135 

9 

999 

7.03 

1.1 

1.0 

13B00 

1000 

28196 

1850 

9999 

9999 

9999 

9999 

9999 

350 

0.15 

0 

13 

4 

125 

0.2 

0 

12 

4 

1530 

0.4 

0 

18 

4 

1830 

1.5 

10 

90193004 

108 

439 

2 

999 

17.2 

1.1 

-1.0 

'14650 

1050 

28196 

1850 

9999 

9999 

9999 

9999 

9999 

(Sheet 
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Table  7  ( 

Contin 

jed)  { 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

1830 

0.5 

0 

1  90193005 

108 

710 

2 

999 

14.8 

1.1 

-1.0 

15825 

1125 

28196 

1850 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

350 

0 

0 

13 

4 

12S 

0 

0 

12 

4 

1830 

1 

10 

90194001 

105 

918 

9 

999 

10.4 

1.6 

1.0 

15800 

1000 

31896 

1850 

9999 

9999 

9999 

9999 

9999 

350 

0.9 

0 

13 

4 

125 

0.35 

0 

12 

4 

1530 

0.5 

0 

18 

4 

1830 

3 

10 

90194004 

SO 

720 

9 

999 

9.09 

1.6 

1.0 

12400 

800 

31896 

1850 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999  . 

9999 

9999 

9999 

9999 

1830 

0 

10 

90195001 

105 

1143 

9 

999 

8.95 

2.6 

1.0 

14850 

850 

35596 

1850 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

1530 

0.5 

0 

18 

4 

1830 

1.5 

10 

90195002 

106 

1141 

9 

999 

8.51 

2.6 

1.0 

19900 

1100 

35596 

1850 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

(Sheet 
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125 

0.45 

0 

12 

4 

1530 

0.56 

0 

18 

4 

1830 

2 

10 

90196003 

140 

1143 

2 

999 

16.6 

2.6 

-1.0 

22950 

1260 

35596 

1850 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

m 

1830 

2 

10 

m 

90195004 

105 

1175 

2 

999 

9.19 

2.6 

1.0 

12100 

700 

35596 

1850 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

9999 

125 

0 

0 

12 

4 

1530 

0 

0 

18 

4 

1830 

0 

10 

Note;  To  convert  to  metric,  multiply  by  these  factors: 
feet  X  0.3048  to  obtain  meters 
feet^  X  0.0929  to  obtain  square  meters 
horsepower  x  9809.5  to  obtain  watts 
inches  x  0.0254  to  obtain  meters 
Description  of  items  in  database: 

First  line: 

1 .  Tost  name 

2.  Vessel  beam,  feet 

3.  Length  of  barges,  feet 

4.  Draft  of  barges,  feet 

5.  Towboat  power,  horsepower 

6.  Vessel  speed  relative  to  ground,  feet/second 

7.  Ambient  velocity  =  Discharge/total  channei  area,  feet/second 

8.  Upbound  =  1.0,  Downbound  -1.0 

9.  Channel  area  left  of  tow  center  line,  feet^ 

10.  Distance  from  left  bank  to  tow  center  line,  feet 

1 1 .  Total  channel  area,  feet' 

12.  Total  water  surface  width,  feet 
Next  1 ,2,  or  4  iines; 

1 .  Distance  from  left  bank  to  velocity  meter,  feet 

2.  Maximum  (upbound)  or  minimum  (downbound)  velocity  alongside  tow,  feeVsecond 

3.  Ambient  velocity  at  this  meter,  feet/second 

4.  Local  depth  at  meter,  feet 

5.  Distance  from  meter  to  bottom,  feet 
Next  1  line: 

1 .  Distance  from  left  bank  to  wave  gauge,  feet 

2.  Water-level  drawdown,  inches 

3.  Local  depth  at  wave  gauge,  feet _ 
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PReOICTED  return  velocitt.  rvsEc 


*****ORL  TESTS  89192-89196 
oqcxxiORL  tests  89198-89203 
TESTS  90172-90179 


OBSERVED  RETURN  VELOCITY.  M/SEC 


Figure  30.  Observed  ORL  prototype  return  velocity  versus  computed  return 
velocity  using  1 995  analytical  method 
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ftoK+fORL  TESTS  89192-89196 
DDDQDORL  tests  89198-89203 
aaaaaORL  tests  90172-90179 
TESTS  90191-90195 


Figure  31.  Observed  ORL  prototype  drawdown  versus  computed  drawdown  using  1995 
analytical  method 
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7  Program  NAVEFF 


The  199S  analytical  method  for  return  velocity  and  drawdown  was 
programmed  in  (^ickBASIC  4.S  language  as  shown  in  the  program  listing  in 
Figure  32.  The  program  prompts  for  metric  or  english  units  and  then  requests 
the  following: 

a.  Total  channel  top  width  in  meters  or  feet. 

b.  Distance  from  tow  center  line  to  left  bank  in  meters  or  feet  (when 
facing  downsteam). 

c.  Total  channel  cross-sectional  area  in  square  meters  or  square  feet. 

d.  Channel  cross-sectional  area  from  tow  center  line  to  left  bank. 

e.  Tow  draft  in  meters  or  feet. 

/.  Total  barge  width  in  meters  or  feet, 

g.  Total  barge  length  in  meters  or  feet.  In  determining  the  displacement 
thickness  used  to  determine  effective  draft  and  beam,  a  single 
temperature  of  17  ”€  is  used  because  of  lack  of  sensitivity  to 
ten:q>erature. 

■h.  Average  channel  velocity  in  meters/sec  or  ft/sec.  The  program  applies  a 
factor  of  1.2  to  the  average  chaimel  velocity  to  determine  the  vessel 
speed  relative  to  the  water. 

i.  Vessel  speed  relative  to  ground  in  meters/second  or  fcet/second. 

j.  Direction  of  travel  U  or  u  for  upbound,  D  or  d  for  downbound. 

The  program  then  outputs  V^(Y)  in  m/sec  or  ft/sec  and  z(Y)  in  meters  or  feet 
at  five  points  on  each  side  of  the  tow.  The  five  points  are  equally  spaced  and 
begin  at  one  tow  width  away  from  the  vessel  center  line  and  end  at  the 
shoreline.  The  program  then  prompts  for  a  new  vessel  speed  using  the  same 
channel  and  tow  geometry.  Tire  program  will  indicate  if  the  speed  entered 
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'PROOMM  HAVBPF.  BAB-BASED  OH  8CHIJF  FOR  BHTIRE  CHANNEL 
'AHD  1995  NEB  RBPT  ON  RETURN  VELOCITY  AND  DRANDONN 
REM  THIS  PRINTS  OUT  ALL  DATA  AMD  6IVBS  A  VISUAL  PICTURE 
CLS 

SCREEN  9 
COLOR  11,  4 

LINE  (700,  600)-(0,  0),  4,  BF 
PRINT 

PRINT  TAB(13};  "PROGRAM  NAVEFF.BAS-SCHIJF  METHOD  PLUS  EMPIRICISM* 

PC  -  0 
PRINT 

PRINT  "  EXPERIMENTAL  DATA  USED  IN  DEVELOPMENT  LIMITED  TO  THE  FOLLOWING:" 
PRINT  *  BLOCKAGE  RATIO  LESS  THAN  B5" 

PRINT  "  DISTANCE  ON  SIDE  OF  VESSEL  GREATER  THAN  10%  OF  TOTAL  CHANNEL  WIDTH" 
PRINT  "  AND  GREATER  THAN  VESSEL  BEAN" 

PRINT  "  VESSEL  LENGTH  GREATER  THAN  40%  OF  CHANNEL  WIDTH" 

PRINT  *  RETURN  VELOCITY  AMD  DRAWDOWN  LIMITED  TO  ONE  BEAM  WIDTH  AWAY" 

PRINT  "  FROM  VESSEL  CENTERLINE  OUT  TO  THE  SHORELINE" 

PRINT  *  VESSEL  SPEED  RELATIVE  TO  WATER  SHOULD  BE  0 . 35-0. 9*LINIT  SPEED" 

FOR  I  >  1  TO  10 

PRINT 

NEXT  I 

PRINT  TAB (25);  "Praia  SPACE  BAR  to  Contlnua" 

SLEEP 

CLS 

t 


FLAG1$  -  "M" 

FLAG3$  -  "N" 

5  INPUT  "ENGLISH  OR  METRIC  UNITS  (E  OR  M)  ",  unl$ 

IF  unl$  -  "E"  OR  uni$  -  "a"  THEN  FLAGl?  -  "Y" 

IF  unl$  -  "M*  OR  uni$  -  «B"  THEN  FLAG2$  -  "Y* 

IF  FLAGIS  •  "Y"  OR  FLAG2$  -  *Y"  THEN  GOTO  7 
GOTO  5 

7  ' 

INPUT  "ENTER  TOTAL  CHANNEL  TOP  WIDTH  ",  BTOTAL 

8  INPUT  "ENTER  DISTANCE  FROM  TOW  CENTERLINE  TO  LEFT  BANK  ",  BLEFT 
INPUT  "ENTER  TOTAL  CHANNEL  AREA  ",  ATOTAL 

INPUT  "ENTER  AREA  LEFT  OF  TON  CENTERLINE  ",  ALEFT 

INPUT  "ENTER  BARGE  DRAFT  ",  D 

INPUT  "ENTER  TOTAL  BARGE  WIDTH  ",  B 

IF  B  <  BLEFT  AND  B  <  BTOTAL  -  BLEFT  THEN  GOTO  9 

PRINT  "DISTANCE  FROM  VESSEL  CENTERLINE  TO  BANK  MUST  BE  >  VESSEL  BEAM" 
GOTO  8 

9  INPUT  "ENTER  TOTAL  LENGTH  OF  BARGES  ",  L 

INPUT  "ENTER  AVERAGE  CHANNEL  VELOCITY,  ^  FOR  UPBOUND,  -  FOR  DOWN  ",  VAN 

10  INPUT  "ENTER  TOW  SPEED  RELATIVE  TO  GROUND  ",  VG 
V  -  VG  +  1.2  •  VAN 

GRAV  •  32.16 
FLAGS  -  "N>' 

IF  uniS  «  "M"  THEN  GRAV  -  9.805 
IF  uni$  -  "■"  THEN  GRAV  -  9.805 
LUNITS  -  "FEET  " 

VUNITS  -  "FEET/SEC* 

AUNITS  -  "SQ  FT" 

IF  unis  -  "M"  OR  uniS  -  "m"  THEN  LUNITS  -  "METERS" 

IF  unis  -  "M"  OR  uniS  -  "n"  THEN  VUNITS  -  "M/SEC  " 

IF  uniS  -  "M"  OR  uniS  "  "m"  THEN  AUNITS  -  "SQ  M  " 


'  SET  WATER  VISCOSITY 


0.0000011  M*«2/SEC  FOR  TEMP  "  17  DEG  C 


Figure  32.  Listing  ^or  PC  program  NAVEFF  (Sheet  1  of  12) 
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VNU  >  .0000011 

IF  LUHIT$  -  "FEET"  THEM  VNU  -  VNU  *  3.2E  *  1.28 

t 

'  CONPUTE  GEOMETRIC  FACTORS 

t 

AH  >  B  •  O 

BRIGHT  -  BTOTAL  -  BLEFT 
ARIGHT  -  ATOTAL  -  ALEFT 
MSIDEL  -  3  *  ALEFT  /  AM 
HSIDER  -  2  *  ARIGHT  /  Kf. 

SKEHL  -  ATOTAL  /  2  /  ALEFT 
BLB  >  BLEFT  /  BTOTAL 

IF  BLB  >  .9  OR  BLB  <  .1  THEN  FLAGS  -  "V" 

SKENR  •  ATOTAL  /  2  /  ARIGHT 
FRB  >  BRIGHT  /  BTOTAL 

IF  BRB  >  .9  OR  LRB  <  .1  THEN  FLAGS  -  "Y” 

IF  FLAGS  >  "N"  THEN  GOTO  20 

PRINT 

PRINT  "  NAPJIIMG,  THE  DATA  ENTERED  FOR  THE  WIDTH  ON  ONE  SIDE  OF  THE  TOW  "; 
PRINT  "  IS  LESS  THAN  10%  OF  THE  TOTAL  CHANNEL  WIDTH.  THIS  IS  OUTSIDE  THE"; 
PRINT  "  LIMITS  OF  THE  EXPERIMENTAL  DATA." 

FOR  I  -  1  TO  IS 

PRINT 

NEXT  I 

PRINT  TAB(2S) ;  "Preas  SPACE  BAR  to  Continua" 

SLEEP 

CLS 

IF  L  >-  .4  "  BTOTAL  THEN  GOTO  14 

CLS 

PRINT 

PRINT  "  WARNING,  THE  DATA  ENTERED  FOR  THE  VESSEL  LENGTH  "; 

PRINT  "  IS  LESS  THAN  40%  OF  THE  TOTAL  CHANNEL  WIDTH.  THIS  IS  OUTSIDE  THE* 

PRINT  "  LIMITS  OF  THE  EXPERIMENTAL  DATA." 

FOR  I  -  1  TO  IS 

PRINT 

NEXT  I 

PRINT  TAB(2S);  "Praas  SPACE  BAS  to  Continue" 

SLEEP 
CLS 
14  ' 

20  H  -  ATOTAL  /  BTOTAL 
N  -  ATO-f  AL  /  AH 
En  -  N 

IF  N  <  85  THEN  GOTO  30 

CLS 

PRINT 

PRINT  "  WARNING,  THE  BLOCKAGE  RATIO  IS  GREATER  THAN  85,  WHICH  IS  OUTSIDE  " 
PRINT  "  THE  LIMITS  OF  THE  EXPERIMENTAL  DATA." 

FOR  I  -  1  TO  15 

PRINT 

NEXT  I 

PRINT  TAB (25) ;  "Prass  SPACE  BAP.  to  Continue" 

SLEEP 
CLS 
30  ' 

REK  'SOLVE  SCHIJF  E&UATION  FOR  DISPLACEMEIIT  CALCULATION 

F 

Z  =  .01 


Figure  32.  (Sheet  2  of  12) 
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40  SCHIJF  -  (1  *  M  *  Z  /  H)  /(N-1-H»Z/H) 

ZT  -  (V  *  2  /  2  /  GRAV)  •  ((SCHIJF  “  2)  +  2  *  SCHIJF) 

Ul  -  V  *  SCHIJF 

IF  ABS((ZT  -  Z)  /  ZT)  <  .00001  THEH  SOTO  SO 
IF  Z  >  20  THEN  GOTO  100 
Z  -  ZT 
GOTO  40 
50  RIM 

REH  COMPUTE  DISPIACEMENT  THICKNESS 
REM 

VDISP  -  V  +  Ul 

RL  -  VDISP  *  L  /  VNU 

DISP  -  .292  *  L  /  (.43429  •  U)G(RL))  *  2. 58 
OE  -  D  4  DISP 
BE  -  B  4-  2  *  DISP 
N  •  ATOTAL  /  BE  /  DE 

REM  'SOLVE  SCHIJF  EQUATION  FOR  RETURN  VELOCITY 

t 

Z  -  .01 

61  SCHIJF  -  (1  +  It  •  Z  /  n)  /  (M-1-N*Z/H) 

ZT  -  (V  “  /  2  /  GRAV)  *  ((SCHIJF  “  2)  +  2  *  SCHIJF) 

Ul  -  V  *  SCHIJF 

IF  ABS((ZT  -  Z)  /  ZT)  <  .00001  THEN  GOTO  52 
IF  Z  >  20  THEN  GOTO  100 
Z  -  ZT 
GOTO  51 
52  REM 

REM  'SOLVE  SCHIJF  EQUATION  FOR  LIMIT  SPEED  USING  MCMTON  RAPHSON 

t 

”LO  -  V 

54  RTEM  •  VLO  2  /  GRAV  /  H 

FV-1-1/N  +  .5*  RTBf  -  1.5  *  RTZM  *  (1  /  3) 

FPV  -  VLO  /  GRAV  /  H  -  (VLO  *  (-1  /  3))  /  (GRAV  *  H)  *  (1  /  3) 

VLN  -  VLO  -  FV  /  FPV 

IF  ABS((VLO  -  VLH)  /  VLO)  <  .0001  THEN  GOTO  55 
VLO  -  VLM 
GOTO  54 

55  VL  »  VLM 

58  VLIMHAT  -  V  /  VL 
PRINT 

IF  VLIMRAT  >  .35  THEN  GOTO  59 

PRINT  •HARNIMO'-*******  SPEED  LESS  THAN  0.35*LIMIT  SPEED  ********  " 

59  IF  VLIMRAT  <  .9  THEN  GOTO  61 

PRINT  "WARNING********  SPEED  GREATER  THAN  0.9*LIMIT  SPEED  *****•*•” 
61  ' 

'  APPLY  CORRECTION  FACTOR 

9 

CF  -  1.7*  -  1,07  *  VLIMRAT 
IF  CF  <  1  THEN  CF  -  11 
Ul  -  CF  *  Ul 

ZT  -  (V  +  Ul)  *  2  /  2  /  GRAV  -  V  *  2  /  2  /  GRAV 
PRINT 

9 

REH  COMPUTE  ■<ALF)  AND  AVERAGE  Vr  FOR  EACH 
REN  SIDE  OF  TOW 
» 


VFACTl 

m 

1.65  -  1.3  •  BLB 

IF  BLB 

> 

.5  THEN  VFACTL  - 

1.35  - 

.7  *  BLB 

VFACTR 

1,65  -  1.3  *  BRB 

IF  BRB 

> 

.5  THEH  VFACTR  - 

1.35  - 

.7  *  BRB 

Figure  32.  (Sheet  3  of  12) 
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VtlAI.  -  U1  *  VFACTL 
VRAR  ••  01  *  VFACTR 
ZSL  -  ZT  *  VFACTL 
ZSR  -  ZT  *  VFACTR 
.ALFL  •  .75  *  N5IDEL  *  .18 
ALFR  -  .75  *  NSIDBR  *  .18 
IF  ALFL  <  1  THEN  ALFL  1 
XF  ALFR  <  1  THEN  ALFR  -  1 
ZALFL  -  ALFL  *  .5 
ZALFR  »  ALFR  *  .5 


VRUI  -  ALFL  •  VRAL 
VRRM  •  ALFR  *  VRAR 
Z8KL  -  ZALFL  *  ZSL 
ZSMR  -  ZALFR  *  ZSR 


PRINT  TAB (22);  "DATA  SHOHN  IS  THE  INPUT  DATA" 

PRINT 

PRINT  "  CHANNEL  TOTAL  AREA  USIIM  V#####.#";  ATOTAL; 

PRINT  •  ■;  AUNIT$; 

PRINT  "  AREA  LEFT  OF  TON  *;  USING  >##/###.#";  ALEFT; 

PRINT  •  AUNIT$ 

PRINT  •  TOTAL  WIDTH  USINO  "tttt/.f;  BTOTAL; 

PRINT  ■  LUHIT$; 

PRINT  •  DISTANCE,  LEFT  BANK  TO  TOW  USING  BLEFT; 

PRINT  ■  LUNIT$ 

PRINT  "  TOW  WIDTH  USING  "##/##.#";  B; 

PRINT  "  LUNIT$;  "  DRAFT  «;  USING  "####.«“;  D; 

PRINT  ■  ■;  LUHIT$ 

PRINT  "  TOW  LENGTH  *;  USING  "#####.#";  L; 

PklNT  "  LDHIT$ 

PRINT  "  TOW  SPEED  RELATIVE  TO  GROUND  USING  •‘/fff.tt"}  VG; 

PRINT  ■  ■;  VONIT$ 

PRINT  "  AVERAGE  CHANNEL  VELOCITY (4-UPBOUND,— DOWN)  USING  •tftf/.ff";  VAN; 

PRINT  ■  VUNIT$ 

PRINT 

t 

FOR  I  -  1  TO  10 

PRINT 

NEXT  I 

PRINT  TAB(25);  "Prass  SPACE  BAR  to  Continua" 

SLEEP 

CLS 

IF  ARIGHT  >  ALEFT  AND  VRRH  >  VRLM  THEN  VRRN  •  VRLM 
IF  ALEFT  >  ARIGHT  AND  VRLN  >  VRRN  THEN  VRLH  -  VRRH 

t 

'CONFUTE  RETURN  VELOCITY  AND  DRAWDOWN  DISTRIBUTION 
/ 

PRINT 

PRINT  TAB(IS);  "CONFUTED  RETURN  VELOCITY  AND  DRAWDOWN  DISTRIBUTION" 

PRINT 

PRINT  TAB(IS) ;  "DISTANCE";  TAB(37);  "RETURN";  TAB(54];  "DRAWDOWN" 

PRINT  TAB(IS);  "FROM  TOW  CL";  TAB(37);  "VELOCITY" 

PRINT  TAB(15);  LUNIT$;  TAB(37);  VUNIT$;  TAB(54) ;  LUN1T$ 

PRINT 

FOR  J  -  1  TO  2 

IF  J  -  2  THEN  GOTO  60 

ALF  >  ALFL;  VRM  -  VRLM 


4 


SALF  •  ZALFL:  ZSH  -  ZSML 
BSIDE  ••  BLEFT 
GOTO  70 

60  ALF  >  ALFR:  VEM  -  VBKH 
ZALF  -  ZALFR:  ZSH  -  ZSMR 
BSIDE  «  BRIGHT 
70  C  -  31  *  LOG(l  /  ALF) 

ZC  -  31  *  LOG(l  /  ZALF) 

f 


■  -BSIDE 

-  B 

/  4  *  (BSIDE  -  B)) 


B)  /  (BSIDE  -  B)) 

-  B)  /  (BSIDE  -  B)) 


-  1  THEM  LOCA$  -  " 

-  5  THEM  LOCA$  «  " 
(J  -  1)  THEN  VIQL 
(J  -  2)  THEN  VYQR 


LEFT  BANK" 
RIGHT  BANK" 
VY!  ZLQ  -  ZY 
ZY 


USING 

TTSING 

USING 


Y; 

"«######.«#";  VY; 
"tttfittm.ttt”;  ZY 


YLAST  -  0 
FOR  I  -  1  TO  5 
IF  J  -  1  THEN  dun 
IF  J  •  2  THEN  dun 
Y  -  (dun  +  (I  -  1) 
i-y  -  ABS(Y) 

VY  -  VRM  •  EXP(C  *  (yy 
ZY  -  ZSH  *  EXP(ZC  •  (yy 
LOCaiS  -  "  " 

IF  J  >  1  AMO  I 
IP  J  -  2  AMD  I 
IF  (I  -  1)  AMD 
IF  (I  -  1)  AMD 
PRINT  LOCAS; 

PRINT  TAB(15); 

PRINT  TAB (30); 

PRINT  TAB (45) ; 

YLAST  -  Y 
NEXT  I 
MEKT  J 

FOR  I  -  1  TO  5 

PRINT 

NEXT  I 

PRINT  TAB(25);  "Press  SPACE  BAR  to  Continue" 

SLEEP 

GOTO  130 

100  PRINT  "SPEED  EQUAL  TO  OR  GREATER  THAN  VLINIT" 
GOTO  10 
130  CLS 
PRINT 
INPUT  ■ 

IF  PRNl  ■ 

PRINT 
INPUT  " 

IF  ANSI 
CLS 
PRINT 
INPUT  ■ 

IF  ANS2 


160 


ENTER  1  FOR  A  PRINTOUT  OF  THE  INPUT  DATA  AMD  RESULTS";  PRNl 
THEN  GOSUB  PAPER 

ENTER  1  FOR  A  VISUAL  DISPLAY  OF  THE  RESULTS";  ANSI 
1  THEN  GOSUB  PICTURE 


ENTER  1  FOR  A  NEW  VESSEL  SPEED  OR  2  TO  QUIT";  ANS2 
THEN  CLS  :  GOTO  10  ELSE  CLS  :  END 


PICTURE: 

CLS 

SCREEN  9 


REM  ADJUSTMENT  FACTOR 
RAT  ■«  560  /  BTOTAL 

OHAX  -  (.75)  •  ((ALEFT  /  BLEFT)  +  (ARIGHT  /  BRIGHT)) 
DNRAT  -  DMAX  *  RAT 
SELECT  CASE  DMRAT 
CASE  IS  >  60 
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■  »  1 

CASE  40  TO  60 

■  >  3 

CASE  26  TO  39 
B  »  3 

CASE  20  TO  2S 

a  >  4 

CASE  1  TO  19 
B  ••  5 
END  SELECT 


RLH  BARGE 
DM  -  D  •  B 
BLRAT  -  BLEFT  *  RAT 
BRRAT  -  BRIGHT  •  RAT 
DRAT  •  D  •  RAT  *  B 
DRTOT  -  4.6 

IF  AUNIT$  -  "SQ  FT*  THEM  DRTOT  -  3.28  *  DRTOT 
DT  -  (DRTOT  -  D) 

DTRAT  -  DT  *  RAT  *  B 

bRAT  -  B  *  RAT 

BLRATA  -  BLRAT  +  40 

blbc  •  BLRATA  -  (bRAT)  /  2 

brbc  -  BLRATA  -t-  (bRAT]  /  2 

brtc  •  brbc 

bite  -  blbc 

blbr  •  160  +  DRAT 

brbr  "  blbr 

brtr  -  160  -  DTRAT 

bltr  “  brtr 

LIME  (bite,  bltr) -(brbc,  brbr),  ,  B 
LINE  (40,  16a)-(600,  160) 


REM  FRAME  WORK 
RAT  -  560  /  BTOTAL 
BLRAT  -  BLEFT  •  RAT 
BRRAT  **  BRIGHT  *  RAT 

IF  VYQL  >  VYQR  THEM  WQ  »  WQL  ELSE  VYQ  -  VVQR 
IF  ZLQ  >  ZRQ  THEN  ZQ  •  ZLQ  ELSE  ZQ  -  ZRQ 
IF  VYQ  <  .2  THEN  VRAX  ■  350:  VP  -  .2 
IF  (VYQ  >  .2)  AND  (VYQ  <  .5)  THEN  VRAT  ■  140:  VP  -  .6 

IF  (VYQ  >  .5)  AND  (VYQ  <  l!)  THEN  VRAT  -  70:  VP  ■  1! 

IF  (VYQ  >  11)  AND  (VYQ  <  2!)  THEN  VRAT  -  35:  VP  -  21 

IF  (VYQ  >21)  AND  (VYQ  <  51)  THEN  VRAT  ■  14:  VP  •  51 

IF  (VYQ  >  51)  AND  (VYQ  <  101)  THEN  VRAT  -  7:  VP  «  10! 
IF  (ZQ  <  .1)  THEN  DRAT  -  700:  DP  >  .1 
IF  (ZQ  >  .1)  AND  (ZQ  <  .2)  THEN  DRAT  -  350:  DP  >  .2 

IF  (ZQ  >  .2)  AND  (ZQ  <  .5)  THBI  DRAT  -  140:  DP  -  .5 

IF  (ZQ  >  .5)  AND  (ZQ  <  1)  THEN  DRAT  >  70:  DP  -  11 

IF  (ZQ  >  11)  AND  (ZQ  <  2)  THEN  DRAT  -  35:  DP  -  21 

IF  (ZQ  >  21)  AND  (ZQ  <  5)  THEN  DRAT  >  14:  DP  51 

PRINT'  VP;  « 

PRINT  TAB(33);  "RETURN  VELOCITY" 

PRINT  TAB (12);  "LEFT";  TAB (60);  "RIGHT" 

PRINT 


IF  uni$  -  "B"  OR  unl$  -  "M"  THEN 
PRINT  "M/SEC" 

PRINT 

PRINT  *  O  «;  BLEFT;  "  <  O’;  TAB(60)  ;  "0  >  ";  BRIGHT 

PRINT  TAB (36) ;  "DRAWDONN" 


4 
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PRINT 
PRINT  "  M" 

ELSE 

PRINT  "F/SEO" 

PRINT 

PRINT  "0  BLEFT;  "  <  0";  TAB (60);  "0  >  ”!  BRIGHT 

PRINT  TAB(36);  "DRAMOOHN" 

PRINT 

PRINT  «  FT" 

END  IF 
PRINT 
PRINT  DP 

LINE  (40,  5)-(40,  160) 

LINE  (40,  80)-(e00,  80) 

LINE  (600,  5) -(600,  160) 

LINE  (40,  290) -(600,  290) 

LINE  (35,  10)-(45,  10) 

LINE  (35.  150)-(45,  150) 

FOR  I  -  1  TO  8 
PRINT 
NEXT  1 

REM  XSECTION 

DMAX  -  (.75)  *  ((ALEFT  /  BLEFT)  +  (ARIGHT  /  BRIGHT)) 

DLEFT  -  ((2  *  ARIGHT)  /  BRIGHT)  -  (DMAX  /  2) 

DRIGHT  -  ((2  *  ARIGHT)  /  BRIGHT)  -  (DMAX  /  2) 

BRRAT  •  BRIGHT  *  RAT 
DLRAT  -  DLEFT  *  RAT 
DPRAT  -  DRIGHT  *  HAT 
DMRAT  -  DMAX  *  RAT 
SELECT  CASE  DMRAT 
CASE  IS  >  60 
B  -  1 

CASE  40  TO  60 
a  ••  2 

CASE  26  TO  39 
a  «  3 

CASE  20  TO  25 

a  -  4 

CASE  1  TO  19 
a  “  5 

END  SELECT 
DLRAT  -  DLRAT  *  B 
DRRAT  -  DRRAT  •  B 
DMRAT  -  DMRAT  *  H 
BLPIC  -  BLRAT  /  2  +  40 
BLP2C  >  BLEFT  *  RAT  +  40 

BRPIC  >  BRRAT  /  2  4-  40  +  BLP2C 

BLPIR  -  DLRAT  *  160 

BLP2R  -  DMRAT  160 

BRFIR  -  DRRAT  *  160 

LINE  (1,  152)-(30,  152) 

LINE  (30,  152)-(40,  160) 

LINE  (40,  160) -(BLPIC,  BLPIR) 

LINE  (BLPIC,  BLPIR) - (BLP2C,  BLP2R) 

LINE  (BLP2C,  BLP2R) - (BRPIC,  BRPIR) 

LINE  (BRPIC,  BRPlR)-(600,  160) 

LINE  (600,  160)-(610,  152) 

LINE  (610,  152)-(640,  152) 

LINE  (BLP2C,  150)-(BLF2C,  5) 
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DLRAT  -  OLEFT  •  BAT 

DRRAT  -  DRIGHT  *  RAT 

DMRAT  -  DMAX  *  RAT 

PRINT  TAB (30);  "BSD  CROSS-SECTION" 

PRINT 

PRINT  "  Bed  Depth  fc  Barge  Draft  are  Increased  by  a  Factor  of";  m; 

PRINT  "for  Display  Purposes" 

PRINT  TAB(lo);  "Barge  Location  and  Size  are  to  Scale(hlt  space  bar  to  continue) 
REM  PROP 

IF  BTOTAL  <  150  THEN 
DIM  PX(15) ,  Py(15) 

IP  PC  -  0  THEN 
FOR  I  >  1  TO  12 
READ  PX(I) 

READ  Py(I) 

NEXT  I 

DATA  10,13,15,17,20,20,24,27,28,30,30,28,28,24,20,20,18,16,15,12,14,10,10,13 
PC  »  1 
END  IF 

FOR  I  -  1  TO  11 

LINE  (PX(I)  +  BLRAT  +  20,  PY(I)  +  140)-(PX(1  +  1)  +  BLRAT  +  20,  Py(I  +  1)  +  14 
NEXT  I 
END  IF 

RD(  DRAHDOHN  t  VELOCITY  LEFT  SIDE 
ALF  >  ALFL:  VRM  -  VRLM 
ZALF  -  ZALFL:  ZSM  -  ZSHL 
BSIDE  -  BLEFT 
C  -  31  *  LOG(l  /  ALF) 

ZC  -  3)  *  LOG(l  /  ZALF) 

r 

'REM  Y-DISTANCE  VY-VELOCITY  ZY-DRAWDOWN 

REM  YL(I)-YLC(I)  VYL(I) -VYLR(I)  ZL(I)-ZLR(I) 

FOR  I  -  1  TO  5 
dun  x  -BSIDE 

Y  >  (dun  4  (I  -  1)  /  4  *  (BSIDE  -  B)) 

YL(I)  -  -Y 

YLC(I)  -  (BLRAT  +  40)  -  (RAT  •  YL(I)) 
yy  »  ABS(YL(I)) 

VYL(I)  ■  VRM  *  EXP(C  *  (yy  -  B)  /  (BSIDE  -  B) ) 

VYLR(I)  -  80  -  (VYL(I)  "  VRAT) 

ZL(I)  X  ZSM  *  EXP(ZC  *  (yy  -  B)  /  (BSIDE  -  B) ) 

ZLR(l)  X  ZL(I)  *  DRAT  4  85 
NEXT  I 

REM  DRAHDONN  t  VELOCITY  RIGHT  SIDE 
ALF  X  aLFR:  VRM  -  VRRM 
ZALF  -  ZALFRi  ZSM  >  ZSMR 
BSIDE  -  BRIGHT 
C  »  3!  *  LOG(l  /  ALF) 

ZC  -  31  *  LOG(l  /  ZALF) 

i 

REM  Y-DISTANCE  VY-VELOCITY  Z¥«DRANDOHN 
REM  YR(I)-YRC(I)  VYH(I) -VYRR(I)  ZR(I)-ZRR(I) 

FOR  I  X  1  TO  5 
dun  X  B 

YR(I)  -  (dun  +  (I  -  1)  /  4  *  (BSIDE  -  B) ) 

YRC(I)  -  YR(I)  *  RAT  4  BLRAT  4  40 
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yy  -  ABS(yR(l)) 

VYR(I)  »  VRM  *  EXP(C  *  (yy  -  B)  /  (BSIDE  -  B) ) 
VYRR(I)  -  80  -  (VYR(I)  *  VRAT) 

ZR(I}  -  ZSK  *  EXP(ZC  *  (yy  -  B)  /  (BSIDE  -  B) ) 
ZRR(I}  -  ZR(I)  *  DRAT  «  85 
NEXT  I 

REM  LOAD  LINES  INTO  LINES.  1  FILE 
OPEN  "LINES.  1"  FOR  OUTPUT  AS  fl 
FOR  I  -  1  TO  5 
IF  I  -  1  THEN 
PRINT  #1.  "BN"; 

PRINT  #1,  USING  "/##";  YLC(I)  ; 

PRINT  j^l  *<  **  • 

PRINT  Fl'  USING  "###";  VYLR(I); 

END  IF 

IF  I  >  1  <  5  THEN 
PRINT  #1,  "  M"; 

PRINT  #1,  USING  "tit”:  YLC(I); 

PRINT  #1, 

PRINT  FI,  USING  "FFF";  VYLR(I); 

END  IF 

IF  I  -  5  THEN 
PRINT  FI  "  H" * 

PRINT  Fl!  USING  "FFF";  YLC(I); 

PRINT  Fi, 

PRINT  FI,  USING  "FFF";  VVLR(I) 

END  IF 
NEXT  I 

FOR  I  -  1  TO  5 
IF  I  -  1  THEN 
PRINT  Fl,  "BM"; 

PRINT  Fl,  USING  "FFF";  YLC(I); 

PRINT  Fl, 

PRINT  Fl,  USING  "FFF";  ZLR(l); 

END  IF 

IF  I  >  1  <  5  THEN 
PRINT  Fl  "K" * 

PRINT  Fl!  USING  "FFF";  YLC(I); 

PRINT  Fl, 

PRINT  Fl,  USING  "FFF";  ZLR(I); 

END  IF 

IF  I  -  5  THEN 
PRINT  Fl,  "M"; 

PRINT  Fl,  USING  "FFF";  YLC(I); 

PRINT  Fl, 

PRINT  Fl,  USING  "FFF";  ZLR(I) 

END  IF 
NEXT  I 

FOR  I  -  1  TO  5 
IF  I  -  1  THEN 
PRINT  Fl,  "BH"; 

PRINT  Fl,  USING  "FFF";  YRC(I); 

PRINT  Fl, 

PRINT  Fl,  USING  "FFF";  VYRH(I); 

END  IF 

IF  I  >  1  <  5  THEN 
PRINT  #1  "M"; 

PRINT  Fl!  USING  "FFF";  YRC(I); 

PRINT  Fl, 
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PRINT 
END  IP 

01, 

USING 

"###"; 

VYRR(I) ; 

IF  I  « 

S  THEN 

PRINT 

#1, 

«M"; 

PRINT 

#1, 

USING 

YRC(I) ; 

PRINT 

#1, 

It  « • 

.  . 

PRINT 
END  IF 

#1, 

USING 

"###"; 

VYRR(I} 

NEXT  I 

FOR  I  - 

1  TO  5 

IF  I  • 

1  THEN 

PRINT 

#1, 

"BM"; 

PRINT 

USING 

"W"; 

YRC(I); 

PRINT 

/I, 

n  H  • 

.  . 

PRINT 
END  IF 

#1. 

USING 

•‘000"; 

ZRR(I)  ; 

IF  1  < 

I  > 

5  THEN 

PRINT 

/I. 

"M"; 

PRINT 

/I. 

USING 

"000"; 

YRC(I)  ; 

PRINT 

01. 

H  M  • 

.  . 

PRINT 
END  IF 

#1, 

USING 

"000"; 

ZHR(I)  ; 

IF  I  - 

S  THEN 

PRINT 

01. 

"M»; 

PRINT 

01. 

USING 

"000"; 

VRC(I); 

PRINT 

01, 

n  H  • 

.  . 

PRINT 
END  IF 

01, 

USING 

"000"; 

ZBR(I) 

NEXT  I 
CLOSE  /! 

1 

a 

1 

LINES  THROUGH  POINTS 

OPEN  "LIMES.!"  FOR  INPUT  AS  /! 
LINE  INPUT  #!,  LVEL$ 

DRAW  "X"  +  VARPTR$(LVEL$) 

LINE  INPUT  /!,  L0RAH$ 

DRAN  "X"  4  VARPTR$(LORAW$) 

LIME  INPUT  #1,  RVEL$ 

DRAW  "X"  4  VAHPTR$(RVEL$) 

LINE  INPUT  #!,  RORAW$ 

DRAW  "X"  4  VAHPTR$(RORAW$) 
CLOSE  /! 

SLEEP 

RETURN 


paper: 

PRINT 

PRINT  "  INSURE  YOUR  PRINTER  IS  ON  WITH  PAPER  IN  IT" 

FOR  I  >  1  TO  4 

LPRINT 

NEXT  I 

LPRINT  TAB(12);  "PROGRAM  NAVEFF.BAS-SCHIJF  METHOD  PLUS  EMPIRICISM" 

LPRINT 

LPRINT  "  EXPERIMENTAL  DATA  USED  IN  DEVELOniENT  LIMITED  TO  THE  FOLLOWING:" 
LPRINT  "  BLOCKAGE  RATIO  LESS  THAN  85* 

LPRINT  "  DISTANCE  ON  SIDE  OF  VESSEL  GREATER  THAN  10%  OF  TOTAL  CHANNEL  WIDTH" 
LPRINT  "  AND  GREATER  THAN  VESSEL  BEAM" 

LPRINT  "  VESSEL  LENGTH  GREATER  THAN  40%  OF  CHANNEL  WIDTH” 

LPRINT  *  RETURN  VELOCITY  AND  DRAWDOWN  LIMITED  TO  ONE  BEAM  WIDTH  AWAY" 
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LPRINT  *  FROM  VESSEL  CENTERLINE  OUT  TO  THE  SHORELINE" 

FOR  I  -  1  TO  8 
LPRINT 
NEXT  I 
LPRINT  T 
LPRINT 
LPRINT  " 

LPRINT  " 

LPRINT  " 

LPRINT  ■ 

LPRINT  ■ 

LPRINT  " 

LPRINT  " 

LPRINT  " 

LPRINT  " 

LPRINT  " 

LPRINT  " 

LPRINT  ■ 

LPRINT  * 

LPRINT  " 

LPRINT  " 

LPRINT  " 

LPRINT  " 

LPRINT 

IF  VLINRAT  >  .35  THEN  GOTO  355 

LPRINT  "WARNING********  SPEED  THRO  WATER  LESS  THAN  0.35*LIMIT  SPEED  ******** 
355  IF  VLINRAT  <  .9  THEN  GOTO  257 

LPRINT  "WARNING********  SPEED  THRU  WATER  GREATER  THAN  0.9*LIHIT  SPEED  ****** 
257  REN 

FOR  I  -  1  TO  5 

LPRINT 

NEXT  I 

LPRINT  TAB(15);  "COMPUTED  RETURN  VELOCITY  AND  DRAWDOWN  DISTRIBUTION" 

LPRINT 

LPRINT  TAB(IS);  "DISTANCE";  TAB{37);  "RETURN";  TAB(54);  "DRAWDOWN" 

LPRINT  TAB (15);  "FROM  TOM  CL";  TAB{37)  ;  "VELOCITY" 

LPRINT  TAB (IS) ;  LUHIT$;  TAB(37);  VUNIT$;  TAB(54);  LUNIT$ 

LPRINT 

FOR  J  -  1  TO  2 
IF  J  >  2  THEN  GOTO  2«0 
ALF  >  ALFL:  VRM  -  VRLH 
ZALF  -  ZALFL;  ZSN  <•  ZSHL 
BSIOE  -  BLEFT 
GOTO  370 

260  ALF  •  ALFR;  VRN  •>  VRRM 
ZALF  >  ZALFR:  ZSN  -  ZSMR 
BSIDE  -  BRIGHT 
270  C  «  31  *  LOG(l  /  ALF) 

ZC  -  31  *  L0G(1  /  ZALF) 

4 


YLAST  •  O 

FOR  I  -  1  TO  5 

IF  J  -  1  THEN  dUB  -  -BSIDE 

IF  J  -  3  THEN  dVlB  -  B 

Y  -  (dUB  +  (I  -  1)  /  4  *  (BSIDE  -  D)) 

yy  -  ABS(Y) 

VY  -  VRM  *  EXP(C  *  (yy  -  B)  /  (BSIDE  -  B)  ) 
ZY  -  ZSM  *  EXP(ZC  *  (yy  -  B)  /  (BSIDE  -  B)) 
LOCA$  -  "  « 


IB(35)  ;  "DATA  SHOWN  IS  THE  INPUT  DATA" 

CHANNEL  TOTAL  AREA  ";  USING  ATOTAL; 

";  AONITZ; 

AREA  LEFT  OF  TOW  «;  USING  "FFFF##.#";  ALEFT; 

";  AUNITS 

TOTAL  WIDTH  ";  USING  BTOTAL; 

" *  LUNIT$  * 

distance!  LEFT  BANK  TO  TOW  ";  USING  •fUff.f"!  BLEFT; 

";  LUNITZ 

TOW  WIDTH  ";  USING  ’‘iitit.fs  B; 

";  LUKITZ;  •  DRAFT  ";  USING  "##«.#/";  D; 

";  LUNITZ 

TOW  LENGTH  ";  USING  "/###/.#";  L; 

";  LUNITZ 

TOW  SPEED  RELATIVE  TO  GROUND  ";  USING  "tffi.tf";  VC; 

";  VUNITZ 

AVERACS  CHANNEL  VELOCITY (-MtUPBOUND, -^DOWN)  »;  USING  "F##F.F#";  VAN; 
";  VUNITZ 
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Figure  32.  (Sheet  12  of  12) 

exceeds  the  Schijf  return  velocity  and  prompts  for  a  new  speed.  At  the  end  of 
each  tabular  output,  the  program  produces  a  screen  image  of  the  cross  section, 
tow,  and  computed  return  velocity  and  drawdown. 

Example  Problem  (based  on  WC  Norman  tow  on  Kampsville  trip  1  field 
data: 

a.  Determine  return  velocity  and  drawdown  distribution  for  a  navigation 
channel  having  a  3S9-m  top  width  and  1,309-sq'm  cross-sectional  area. 
The  tow  is  222  m  from  the  left  bank  which  results  in  a  cross-sectional 
area  left  of  the  tow  of  8(X)  sq  m.  The  barges  draft  2.74  m,  have  a  total 
beam  of  32  m,  and  a  total  length  of  barges  of  238  m.  The  downbound 
tow  travels  at  2.9  m/sec  relative  to  ground  against  an  average  ambient 
velocity  of  0.49  m/sec. 

b.  Solution:  Enter  an  "m"  to  indicate  the  use  of  metric  units  followed  by 
the  geometric  factors.  Enter  an  average  channel  velocity  of 

-0.49  m/sec,  since  this  is  a  downbound  tow.  Enter  a  vessel  speed 
relative  to  ground  of  2.9  m/sec.  The  vessel  length  is  used  to  ensure  the 
ratio  of  vessel  length/chatmel  width  exceeds  0.4  and  to  compute  the 
displacement  thickness  to  determine  the  effective  draft  and  beam  of  the 
vessel.  The  program  output  is  shown  in  Figure  33. 
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PROGRAM  HAVKFF. BAS-SCHIJF  METHOD  PLUS  EMPIRICISM 


EXPERIMENTAL  DATA  USED  IN  DEVELOPMENT  LIMITED  TO  THE  FOLLOHIHG; 
BLOCKAGE  RATIO  LESS  THAN  8S 

DISTANCE  ON  SIDE  OF  VESSEL  GREATER  THAN  10«  OF  TOTAL  CHANNEL  WIDTH 
VESSEL  LENGTH  GREATER  THAN  40%  OF  CHANNEL  WIDTH 
RETURN  VBLOCITV  AND  DRAWDOWN  LIMITED  TO  ONE  BEAM  WIDTH  AWAY 
FROM  VESSEL  CENTERLINE  OUT  TO  THE  SHORELINE 


DATA  SHOWN  IS  THE  INPUT  DATA 

CHANNEL  TOTAL  AREA  1309.0  SQ  M  AREA  LEFT  OF  TOH  800.0  SQ  H 

TOTAL  WIDTH  359.0  METERS  DISTANCE,  LEFT  BANK  TO  TOW  222.0  METERS 

TOH  WIDTH  32.0  METERS  DRAFT  2.74  METERS 

TOH  LENGTH  238.0  METERS 

TOW  SPEED  RELATIVE  TO  GROUND  2.90  H/SBC 

AVERAGE  CHANNEL  VELOCITY (T>UPBOUKD,->DOHN]  -0.49  M/SEC 


COMPUTED  RETURN  VELOCITY  AND  DRAWDOHV  DISTRIBUTION 


DISTANCE 

RETURN 

DRAWDOWN 

PROM  TOH  CL 

VELOCITY 

METERS 

M/SEC 

METERS 

LEFT  BANK  -222.0 

0.155 

0.049 

-174.  sS 

0.185 

0.053 

-127.0 

0.220 

0.058 

-79.5 

0.263 

0.064 

-32.0 

0.313 

0.070 

32.0 

0.364 

0.084 

58.3 

0.324 

0.079 

84.5 

0.289 

0.075 

110.8 

0.257 

0.071 

RIGHT  BANK  137.0 

0.229 

0.067 

Figure  33.  Output  from  PC  program  NAVEFF 


Chapter  7  Program  NAVEFF 


85 


Results  and  Conclusions 


8 


The  recommended  1995  analytical  method  presented  herein  is  based  on 
conservation  of  energy  plus  empiricism  to  define  the  distribution  of  return 
velocity  and  drawdown.  The  1995  analytical  method  is  summarized  to 
compute  the  following: 

a.  Schijf  return  velocity  and  drawdown  using  Equations  1  and  2. 

b.  Average  return  velocity  using  Equation  13. 

c.  Compute  average  return  velocity  on  each  side  of  the  ves-rel  using 
Equations  14  and  15. 

d.  V„/V„  “  a’  asing  Equation  16. 

t.  C  using  Equation  17. 

/.  V,(Y)  using  Equation  10.  V,(Y)  is  linearly  added  to  or  from  ambient 
currenis. 

g.  z(Y)  using  Equations  18  through  22. 

The  1995  analytical  method  should  be  limited  to: 

a.  N  l«ss  than  85.  While  the  vast  majority  of  the  data  had  N<  52,  the 
limited  data  from  the  Clark’s  Ferry  pool  572.7  physical  model  tests  for 
N  of  about  85  resulted  in  significant  scatter  but  exhibited  the  correct 
trend  about  the  line  of  perfect  agreement. 

b.  Tow  length  greater  than  40  percent  of  the  channel  width. 

c.  Distance  on  both  sides  of  tow  center  line  equal  to  or  greater  than 
10  percent  of  the  total  channel  width. 

d.  Distance  Y  from  the  tow  center  line  greater  than  tow  width  B. 

e.  Vessel  qieed  equal  to  35  to  90  percent  of  the  limiting  speed. 


Chaptar  8  Raiultt  and  Condutlont 


/  The  predictive  method  presented  herein,  is  applicable  to  river  reaches 
that  can  be  characterized  by  a  single  cross  section.  One  would  not 
expect  these  techniques  to  provide  valid  results  at  the  end  of  an  island 
or  in  other  areas  where  the  cross  section  varies  rapidly. 

Areas  of  needed  research  for  this  method  are  more  data  for  N  greater  than 
52,  variable  tow  length  less  than  0.4  times  the  channel  width,  and  better  data 
supporting  the  distribution  of  water-level  drawdown.  Future  versions  of  the 
NAVEFF  model  will  include  propeller  jet  ejects  as  well  as  short-period  wave 
activity. 
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Appendix  A 
Notation 


Undisturbed  channel  area 

A^  Disturbed  channel  area  around  midsection  of  vessel 

A^  Submerged  cross-sectional  area  of  vessel  at  midsection 

0  Coefficient  that  varies  with  V/V^  in  the  Schijf  equations  to  improve 

agreement  between  observed  and  computed  return  velocity  and 
drawdown. 

£«  Chaimel  width 

B  Beam  of  vessel 

Distance  from  vessel  center  line  to  shoreline 
g  Gravitational  acceleration 

h  Undisturbed  average  water  depth 

L  Vessel  length 

N  Ratio  of  channel  cross-sectional  area  to  vessel  cross-sectional  area 
V  the  vessel  speed 

Vi  Maximum  speed  (called  the  critical  or  limiting  speed) 

V,  Average  return  velocity 

Maximum  return  velocity 

V„  Average  return  velocity  for  each  side  of  vessel 
V^  Maximum  return  velocity  on  each  side  of  the  tow 


Appendix  A:  Notation 


V,(Y)  Maximum  return  velocity  during  tow  event 
Y  Distance  from  center  line  of  vessel 

z  Average  water-level  drawdown 

Z„  Average  drawdown  for  each  side  of  vessel 
Maximum  drawdown  for  each  side  of  vessel 
Z(Y)  Maximum  water-level  drawdown  during  tow  event 


Appendix  A:  Notation 
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